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ABSTRACT 


An activation-contraction model for skeletal muscle is 
developed based on the sliding filament theory and calcium binding. 
The effects on the isometric twitch and tetanus are studied by 
changing the model parameters which include rates of making (f) and 
breaking (g) of cross-bridges, amount of calcium released by a 
Stimulus (Cg), rates of calcium binding (a) and re-uptake (8), and 
series stiffness (k). The characteristics of the isometric twitch 
are most sensitive to changes in f or B. 

In a muscle with a poorly-developed sarcotubular system, 
both Co and B may be reduced. This gives a twitch which has a smaller 
peak tension (Pt), a lower maximum rate of tension development (dP/dt), 
a longer contraction time (Te) and half-relaxation time (Tip). A 
tetanus which is produced by reducing Cg and B has a smaller maximum 
tension which is better maintained, a longer time to peak and a lower 
fusion frequency. Including a threshold in the activation mechanism 
intensifies the above effects and gives a sigmoid frequency-tension 
curve. These effects are observed experimentally. Reducing k gives 
a twitch with a smaller Pz and dP/dt, and a longer To. 

The magnitudes of these effects are compared with experi- 
mental data. The effects of temperature and of a tetanus on the mech- 
anical properties of muscle are not well enough Known to have been 
included in the model. Experiments have been carried out to study 
these effects more quantitatively using the frog sartorius muscle. 
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decreases by 49% and 22% respectively with an increase in temperature 
from 33° to 30°C. A tetanus of a few seconds duration produces a 
delayed increase in Py when temperature is below 25°C. 

Following a tetanus, Py rises nearly lineariy to a maximum 
and decays with an approximately exponential time course. This is in 
agreement withtrehyoothesis that this phenomenon is The result of a 
fast-decaying inhibition and a slow-decaying potentiation. These 
processes have @19's of 2.7 and 2.9. The slow process also has a half- 
time of about 13 minutes at 20°C. It is suggested that inhibition is 
related to Catt depletion from the release site while potentiation is 
associated with a decreased 8. 

The frequency response of the isometric sartorius muscie 
was determined by analysis of the tension generated during short 
periods of random stimulation of the muscle nerve. It resembles 
closely the Fourier transform of the isometric twitch, and both are 
well represented by the frequency response of a linear second-order 
system. The second-order parameters: zero-frequency gain (Gq), 
natural frequency (f,,) and damping ratio (¢) provide useful descriptions 
of the changes in the muscle response caused by temperature changes or 
tetanus. For example, f, is found to vary linearly with temperature in 
the range 3 to 31°C. 7 is increased following a tetanus. 

At all temperatures, Te, TLp and the ratio Tiy/Tg are 
increased by a tetanus; but the increase in the ratio becomes marked 
at higher temperatures. Analysis using a linear second-order system 


shows that ¢ depends on the ratio Tir/Te rather than the absoiute value 
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CHAPTER 1 


INTRODUCTION. 

Following the introduction of the concepts of active state 
and degree of activation by A.V. Hill (1949), it has been widely 
accepted that a muscle must be activated before it is capable of 
shortening or developing tension. The capacities to shorten and to 
develop tension in turn define the active state of the contractile 
component (Hill, 1951b). It is also believed that force generation 
in a living muscle is due to the continual making and breaking of 
cross-bridges between the overlapping thick and thin filaments 
(Huxley and Niedergerke, 1954; Huxley and Hanson, 1954). A.F. Huxley 
(1957) has worked out a detailed mathematica! formulation of this 
general nature, and has shown that it can account for many features 
of contraction. In his scheme, the making of cross-bridges was 
assumed to be the main rate-limiting factor in steady shortening, 
while breaking after the performance of work by a cross-bridge was 
relatively rapid. The results of this mathematical analysis agree 
well with the 'characteristic equation' formulated by A.V. Hill 
(1938); the latter equation relates speed and load in isotonic 
shortening, and was derived on the basis of muscle energetics. 
However, the results did not fit the transient responses, such as 
the change in tension following a sudden alteration of length or 
vtee versa (Podolsky, Nolan and Zaveler, 1969; Huxley and Simmons, 
Kobra pie Podolsky, Nolan and Zaveler (1969) reversed the assumptions 


about the rate constants, attachment of cross-bridges rapid and 
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breakage rate limiting. In this way, they were able to fit many 
aspects of transient responses but not the thermal data. Recently, 
Huxley and Simmons (1971) applied the theories of visco-elastic 
behavior of high molecular weight polymers to the attachment and 
detachment of cross-bridges and showed that their proposal! could 
account for the successful aspects of both of these theories. 

The basis of the sliding filament model can be used to 
explain many phenomena in muscle and its properties. For example, 
the observation of Gordon, Huxley and Julian (1966) that the maximum 
isometric tension was related to the sarcomere length was explained 
if the number of interacting cross-bridges was proportiona! to the 
amount of overlapping of the filaments. The elasticity which 
accompanied a very small stretch in length in a passive muscle was 
attributed to some 'active' interaction between the filaments (Hil!, 
1968). 

A number of hypotheses were proposed to explain the 
mechanism of force generation at the molecular level. Features 
such as ATP-splitting, calcium ions, electrostatic phenomena, myosin 
and actin molecule structures were usually incorporated into the 
hypotheses. Davies (1968) suggested that each cross-bridge shortened 
by folding at a number of points. Huxley (1969) proposed that the 
head of the myosin molecule rotated relative to the thin filament, 
acting as a lever which pulled the thick filament along by a link 
which was also part of the myosin molecule. The flexible linkage in 
the latter scheme enables equal force production over a wide range 


of interfilament separations. Thorson and White (1969) suggested a 
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hypothesis which related the attachment-detachment cycle of cross- 
bridges to the measured length-tension dynamics of the active insect 
fibrillar flight muscle, in which length change was small. 

The sliding filament model and Hill's 'characteristic 
equation’ are adequate if one requires only the steady-state relation- 
ship, for example, when a muscle is maximally activated. To account 
for the time course of a twitch or tetanus, an activation mechanism 
must be included. Considerable evidence indicates that the binding 
of calcium (Catt) by the myofilaments regulates the process of con- 
traction and relaxation (Sandow, 1965; Weber, 1966; Ebashi, Ebasni 
and Kodama, 1967). Based on this concept, Julian (1969) extended 
the sliding filament model by including an activation mechanism 
which depended on a number of factors. The latter includes the rates 
at which calcium is bound and released, and the amount of calcium 
available. The time course of an isometric twitch can then be 
computed. However, to compute a tetanus or partially-fused twitches, 
summation of effects must be described. Together with the parameters 
in the original sliding filament model, there are no less than eight 
parameters, each of which can be changed independently of the others. 
The existence of this large number of parameters poses some interesting 
problems: 

(1) What are the effects of changes in these parameters on the 
characteristics of an isometric twitch or tetanus? 

(2) How closely can these changes be related to the structural 
differences observed in skeletal muscles? and finally, 


(3) Can these changes be used to explain changes in the char- 
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acteristics of an isometric twitch under different physiological 
conditions, for example, after a tetanus? 

In the event of excitation-contraction coupling in skeletal 
muscle, it is well known that potential change in the plasma membrane 
eventually leads to the release of Catt. Extracellular Catt is also 
required for generation of the muscle action potential (Ishiko and 
Sato, 1957; Edman and Grieve, 1964). Frank (1958, 1960, 1962) found 
that by soaking the toe muscle of the frog in a Catt-free solution 
for a few minutes, he could completely eliminate the potassium- 
induced contracture without reducing the potassium depolarization. 
The rate at which this inhibition developed was determined by the 
rate at which Catt left the extracellular space of the muscle. The 
contractile mechanism of the cells was, however, unaffected by this 
treatment because their response to an appropriate stimulus was not 
changed by the removal! of extracel lular Cas, WTS! eden to 
believe that a 'loosely bound fraction! of Catt was necessary to 
elicit a contracture in response to depolarization of the membrane, 
but another source of Ca*t, the "firmly bound fraction', appeared 
to be needed for contraction itself. The exact relationship between 
these two Catt fractions remains unknown. Discovery of the trans- 
verse tubules (T-tubules) with their close apposition to the sarco- 
plasmic reticulum (Porter, 1956, 1961; Porter and Palade, 1957; 
Bennett and Porter, 1953; Bennett, 1955; Anderson-Cedergren, 1959), 
the calcium accumulation property of the latter structure (Marsh, 
1952; Ebashi, 1961a; Hasselbach, 1964), the continuity of the T- 


tubules with the plasma membrane (Franzini-Armstrong and Porter, 
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1964; Huxley, 1964; Page, 1964; Birks and Davey, 1969) and their 
electrical coupling (Falk and Fatt, 1964; Peachey, 1965), the local 
stimulation experiments of Huxley and Taylor (1958) and Podolsky 
and Costantin (1964), suggest that in the excitation-contraction 
process, potential change of the plasma membrane produces similar 
change in the T-tubules, and probably also in the sarcoplasmic 
reticulum from which Catt is then released. This hypothesis is 
also supported by the observation of Costantin and Podolsky (1966) 
and Lee, Ladinsky, Choi and Kasuya (1966) that local depolarization 
of the triad sacs appears to activate contraction. Activation via 
the T-tubules is probably more important in the fast twitch muscle, 
which has a well-developed internal membrane system compared to the 
Slow muscle (Page, 1968; Flitney, 1971). This is also the 'fast! 
process suggested by Hill (1949) as opposed to the 'slow' diffusion 
process of activation introduced by Heilbrunn and Wiercynski (1947). 
In any case, the final step in activation is an increase in concen- 
tration of calcium in the myofilament space of the cell. 

The extensive distribution of the internal membrane system 
within the cell, and the fact that at least part of this system can 
accumulate calcium, raise the possibility that activation and relax- 
ation in vivo involves intracellular translocation of calcium rather 
than movement across the cell membrane. Direct evidence comes from 
the studies of Winegrad (1965a, b, 1968, 1970). By exchanging some 
of the intracellular Catt with Ca*®, and using autoradiographic 
techniques, he was able to demonstrate the location of Cato inva 


resting and in an active muscle. For example, Ca*> was found to be 
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located mainly in (a) the terminal cisternae in a resting muscle, 
(b) the space between the myofilaments at the peak of a tetanus, 
and (c) in the A-band between the myofilaments (presumably in the 
intermediate cisternae and in the longitudinal tubules) in a 
declining tetanus. The fast exchange between Catt in the terminal 
cisternae and Cat? jn the bathing solution again demonstrated the 
accessibility of the terminal cisternae from the extracel !ular 
space. More recently, Connolly, Gough and Winegrad (1971) studied 
the effect of tetanus on the characteristics of the isometric twitch, 
and found maximal potentiation did not occur until about 10 to 12 
seconds after the tetanus at 20°C. The shape of the isometric 
twitch as characterized by peak tension, maximum rate of tension 
development and contraction time, could be expiained if two processes 
with half-times of about 2 seconds and 2 to 3 minutes occurred after 
a tetanus. The faster process also had a time course similar to that 
of the post-tetanic calcium movement in the sarcoplasmic reticulum. 
The relaxation phase of the twitch and the slow process were, 
however, not studied to any great detail. It seems certain that the 
characteristics of a twitch are governed mainly by the calcium dis- 
tribution within the muscle cell, although a correlation between 
them is difficult. Studies of muscle activation along these iines 
are, however, useful to obtain indirect information on the under- 
lying processes. 

There are a number of other factors which have a direct or 
indirect influence on the characteristics of a twitch. Temperature 


variation is probably the most significant environmental factor 
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which causes changes in the rates of biological reactions. This is 
especially true among poikilothermous animals. The effect of temp- 
erature variation on biochemical and mechanical reactions is complex. 
It would be of value if direct experimental measurements of temperature could 
be compared with the effects of parameter changes in a model set 

up to represent the activation and contraction processes. Hil! 
(1951c) reported the peak twitch tension in frog muscle was decreased 
when temperature was aes in the range 0 to 20°C while tetanus was 
somewhat increased. The fact that the maximum tension in a twitch 
was considerably less than in a tetanus had been attributed to on- 
coming relaxation allowing insufficient time for internal shortening 
to be completed. This would require that the temperature coefficient 
of the velocity of shortening would be substanttally less than that 
of the decay of activity. Earlier, Hartree and Hill (1921) also 
found the temperature coefficient of the maximum rate of fall of 
tension in a twitch was about 1.5 times that of the maximum rate of 
rise of tension. In the neuromuscular system, however, the effects 
of temperature change are usually manifested in the rate change of 
some mechanical processes. 

Hutchinson, Koles and Smith (1970) have reported a linear 
relationship between conduction velocity and temperature in single 
myelinated nerve fibers in Xenopus laevis (the African toad). 

Similar studies on the speed of response of muscle as characterized 
by the time course of twitches due to temperature variation have been 
scarce. Close (1965) and Close and Hoh (1968a) suggested that the 


ditferences in isometric twitch contraction times of different 
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muscles (fast and slow) might be due largely to differences in the 
twitch/tetanus ratio or peak twitch tension. A decrease in temp- 
erature decreased the intrinsic speed of shortening and the rate 

of removal of activator in all these muscles, thereby causing an 
inversely proportional increase in contraction time with little or 
no change in peak tension. In addition, it also increased peak 
tension of fast muscles in much the same way that repetitive stim- 
ulation increased it, by increasing the degree of activation of 
individual muscle fibers with little or no change in the time course 
of the response. Since the processes of activation, contraction 
and relaxation primarily involve the movement and distribution of 
intracellular Catt, it is important to relate them to the mechanical 
response of the muscle, and to temperature. 

The characteristics of an isometric twitch can be studied 
using a model or by experimental investigation. The method of 
Julian (1969) for an activation-contraction model is extended in 
Chapters 2 and 3. The activation mechanism is extended to model 
repetitive stimulation and to include a threshold for contraction. 
Chapter 4 shows the effect of parameter changes of the model on the 
characteristics of the isometric twitch and the tetanus, and the 
rationale of changing any parameter is discussed. In later chapters, 
some effects of tetanus and temperature on the isometric twitch of 
the sartorius muscle of the bullfrog are demonstrated. Results 
obtained from time domain and frequency domain analyses are compared. 
The time course of the mechanical response is related to the speed 
of response of the contractile component and to rates of reactions 


of other processes. 
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CHAPTER 2 


DESCRIPTION AND ASSUMPTIONS OF THE MODEL. 
The components of the model and the major events in an 


activation-contraction cycle are shown in Figure 2.1. 


2.1 Force Generator. 

The basic units of the force generator are essential ly 
those proposed by Huxley (1957) and consist of cross-bridges on 
the myosin filament and specific sites on the actin filament. A 
cross-bridge on the myosin filament is able to oscillate about its 
equilibrium position and to hook up with specific sites on the 
actin filament. Cross-bridges can only be made on one side of the 
equilibrium position. This results in net generation of force in 
one direction and a pulling action on the actin filament. Huxley 
(1957, p. 284) described the . rate of change of interacting 


cross-bridges by a first-order equation: 


a nr) 7 ey FIG Cans 


where f and g are position-dependent rate constants of making and 
breaking of cross-bridges respectively, and m the proportion of 
interacting cross-bridges. Figure 2.2 is redrawn from Huxley 
(1957, p. 282) and shows the magnitude of f and g as a function of 
x, the distance of the thin filament site from the equilibrium 


position of the myosin cross-bridge; h is the positive distance 
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Figure 2.1 A block diagram of the model. Arrows indicate direction 


in which events are assumed to flow. The dashed line 
indicates possible motion feedback path. The circle 
labelled "force" shows the position of the force trans- 


ducer. 
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Figure 2.2 


Dependence of f and g on x (from Huxley, 1957). Dashes 


for f and continuous line for g. 


f1 and g, are res- 


pective values of f andgatax=h. go is value of g 


when zx < 0. 


away from the equilibrium position beyond which a cross-bridge may 
not form a new interaction with a thin filament site and is assumed 
to have a value of 100 A (Podolsky, Nolan and Zaveler, 1969). In 
the resting state, this interaction is prevented by some inhibition 
mechanisms. 

Solving equation (2.1) with the constraint that the con- 
tractile unit is held isometric, gives: 

Seep pete G27 

n(u) ear : eae eek 
Whererumis thegdistance x normalized with respect toh; mn, Is the 
initial n value at any x. As t becomes large, n = f/(f +g), the 
isometric steady-state value for n. 

The force, F, produced by an’ interacting cross-bridge at 
a distance x away from the equilibrium position is assumed to be 
kx, where k is a stiffness coefficient. The force produced by all 


the interacting cross-bridges is then: 


The infinite upper and lower limits on the integral can be taken to 
mean that the integration is extended from the zero position until 
limits are reached at which the value for the integral no longer 
changes appreciably. 

Reducing equation (2.3) to a normalized form suitable for 


the isometric steady-state by substituting = f/(f + g) gives: 
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where Fy is the maximum isometric force, and Kg is a normalizing 
factor for the stiffness coefficient. Since F/Fp = 1 in the iso- 


metric steady-state, k/kg = 2(f + g)/f. Equation (2.3) can now be 


written as: 
ua nudu GA75> 


where P = F/Fy. The equivalent, instantaneous stiffness of the force 


generator is given by: 


kg = See fondu (2.6) 


2.2 > Activation} 
(t) Isometrte Twitch. 

An activation mechanism is necessary when non-steady-state 
processes, such as the time course of tension change in an isometric 
twitch is to be calculated. Following Julian (1969), the rate con- 
stant f which governs the rate at which cross-bridges interact is 
assumed to depend on the binding of calcium. Huxiey's (1957) con- 
straints on f as shown in Figure 2.2 are that f = 0 when u < O or 
u> 1. The value of f is weighed by an activation factor y(t) such 
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y(t) varies between O and 1. 

Experimental studies using the frog sartorius muscle (Jewel | 
and Wilkie, 1958, 1960; Close, 1962) suggest that in an isometric 
twitch produced by a single stimulus, the active state rises to a 
peak in about 30 to 60 msec and decreases to a very tow level by 
about 1.2 sec at 0°C. Jewell and Wilkie (1958) also recorded a 
twitch/tetanus ratio of up to .92, indicating that activation might 
be complete at this temperature. For these reasons, Julian (1969) 
assumed that in a twitch following a stimulus, y rose from zero to a 
maximum 40 msec later, and then decayed back to very near zero 1.2 
sec after the stimulus. 

The time course of activation can be represented by two 
consecutive first-order reactions: 


e] B 


Ce A ry 


where C = precursor; 

A = aetivator: 

E = end-product; 

a = rate constant of binding of calcium by the myofi laments; 

8 = rate constant of release of calcium by the myofilaments and 

re-uptake by the sarcoplasmic reticulum. 

C is assumed to be the amount of calcium released by a stimulus and 
E is a form of Catt which is ready for re-uptake into the sarcoplasmic 
reticulum. Reverse reactions are not allowed. Activation can only 
take effect when A is present. 


The rate equations for these reactions are: 
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al) = -a[C] (2.8) 


24) = atc] - tA] (2.9) 


alt) = BIA] : (2.10) 


where [C] is the concentration of C. 
Assuming that at ¢ = 0, [C] = Co, [A] = [Z] = 0, equations 


(2.9) and (2.10) can be integrated to give the concentration of A: 


Co o 


[4] = 5 


[ exp(-at) - exp(-8t) ] CZet) 


[A], the concentration of the activator, is assumed to be the amount 
of calcium bound to the myofilaments. The following type of function 


is selected to represent y: 


Coa 
rr ar [ exp(-at) - exp(-Bt) ] KZ eal 2 


where [A], = amount of calcium bound at peak of activation; 


ti 


Co amount of calcium released by one stimulus; 

t = time in seconds. 

Julian (1969) assumed that Co was the total calcium in 1 gm 
of muscle and it ws given a value of 2 umoles (Gilbert and Fenn, 1957). 
Using a value of [A]p = .15 umole (Weber, 1966) and 40 msec for the 


time-to~peak activation, the rate constants a and B could be calculated 


from equation (2.12) and its first time derivative. In the present 


. 
i oh Ya sr} 4 ni @ oT ye 7 mrt sf) ; 
nid ‘ ie ; L% : i ‘ rer Rs |b ay , yO : a ioies to 
' , ¢ i saaqva pe er oofoelae € 
i-—_> # A | 
shit ecw { Lee ) 2 - (a4 rf ' — eee Gai, 7 © ¢ 
- - J 3/7 | g . 
a sm oe ES ees ee et Ae iat fe UR Ie TG A AB. ede 
7 : i >» BAX Vi ines a! 2 ae 4D tee =- ink, 4 
aoe te eas - danse oi amit ts Jie 
s bal piinlés Hooks . am o\i tec) tequcee ma } on the 
puna mem. et oe 
nn Be a 
a si) ane: " cia oe 9 y 
ze, = | — —_ apaliis 9° ge outne ant 


: ~~ a Ber > the — 6 
: -" <= on ii . ais od ie i. 33 “apo goT =a 
mr 


° eta ; eoft ony bak (816) naiteu 
hare ws ing his 7 sth: 


16 


model it is assumed that Cg is the amount of calcium released by a 
Single stimulus, instead of the total calcium in 1 gm of muscle. 

a also has a numerical value greater than 8. These assumptions are 
more realistic because only a fraction of the total calcium is used 
in each twitch and binding of calcium by the myofilaments is a 

faster process (Huxley and Simmons, 1971). To keep the time course 
and magnitude of y identical, the values of a and 8 were switched 

and Cg was set to a fraction equal to .192/2 of the total calcium. 
The standard values for the variables governing the activation factor 
are thus [Alp = .15 umole, Co = .192 umole, a = 65/sec, B = 6/sec. 
The time course of the activation factor calculated using these values 


is shown in Figure 2.3. 


(tit) Repetitive Stimulation. 

To account for the time course of the activation factor in 
a fused tetanus, an arbitrary function can be written so that it 
increases from zero to one in a way similar to that given by equation 
(2.12) but remains at one thereafter (Julian, 1969). If the 
frequency of stimulation is low so that an unfused tetanus results, 
the time course of the activation factor fo!lowing each stimulus must 
be calculated. If Ag is the amount of residual activator substance 
as a result of previous activity, then solving equations (2.9) and 
(2.10) gives: 


Co yO Ag 


[ exp(-at) - exp(-8t) ] +> exp(-8t) (2.13) 
[Alp 
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Figure 2.3 
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Postulated time course of the activation factor, y, 
for isometric twitch (equation (2.12) ). 
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where Co ln is the amount of calcium released by the wt stimulus. 


cr Sy = amount stored before the wth stimulus, 


a, woe 


Sy = Sy-1 - Co,n=1 


Solving these equations gives: 


x 92 te ¥ POP mys 
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and: Com = .192 (1 - 28 yt (2.14) 
1 


From equation (2.14), it is obvious that Co,w is greatest 
in the first stimulus. Any subsequent values are decreased as the 
number of stimuli is increased. The amount of calcium that remains 
unbound immediately before a stimulus (other than the first) is smal | 
and is neglected in the computation. This is because any free Cae 
is readily bound (with a rate a = 65/sec), as can be seen from 
equation (2.8). 

The inclusion of the factor Ag in equation (2.13) enables 
y to have values greater than 1 under certain circumstances. However, 


the tension produced can only increase to a maximum of about 1.2 for 


Ag 
reasons which will be given later. In fact Faz has a value equal to 
p 
y immediately before the next stimulus is delivered. Its magnitude 


depends on the previous value of Ag and the time lapse since the 
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previous stimulus. 


(ttt) Threshold. 

The necessity to include a threshold for developing tension 
follows from the fact that in slow muscle fibers, a single stimulus 
does not elicit a twitch; considerable tension can only be generated 
by repetitive stimulation (Kuffler, Laporte and Ransmeier, 1947; 
Kuffler and Vaughan-Williams, 1953). Studies on the mechanical 
threshold of siow muscle are very few. Although a strict comparison 
is impossible, it is important to have an approximate value of Catt 
concentration ([Catt]) for the threshold. Weber and Winicur (1961) 
found a threshold [Cao] of about .2 umole for a 'natural! actomyosin 
preparation and Portzeh!l, Caldwell and Ruegg (1964) reported a 
value of .3 to 1.5 umole for muscle fibers of the crab. 

Julian (1971) studied the effect of calcium on the force- 
velocity relation of briefly glycerinated frog muscle fibers. He 
observed the contractions induced and the steady force developed at 
about 4°C by lowering the pCa in the relaxing solutions to various 
levels determined by the ratio of calcium and EGTA added. The 
relationship between steady force generated and pCa was S-shaped and 
very steep. Contractions were never observed above pCa 7. The data 
for his Force-pCa curve (Julfan, 1971, p. 130) are replotted in 
Figure 2.4 with [Can] instead of pCa. Below a Casal of .12 umole, 
no force is produced; and the steady force rises from zero to half its 
maximum value with a change of Teer of .05 umole. 


To find the value of y which corresponds to [eae] = .12 umole, 
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Figure 2.4 
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Tension - (Gae] curve for determination of threshold 
[Catt] for contraction. (Data from Julian, 1971). 
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the value of in equation (2.2) is considered. At the isometric 
steady-state, n = fy/(f, + g,) = .8125, using values of f) = 
81.25/sec and g, = 18.75/sec (Julian, 1969). The steady force P 


is related to y as follows: 


Suppose f = 81.25 y (so that when y = 1, f = f) = 81.25) 


23 5p 


In Figure 4; when P = 3, n = f/(f +91) = + or f = 91, assuming the 


steady-state value of m at any P depends only on f. Substituting 


This, value, for .into equation .( 2.15): 


(2.16) 


where Yip is the value of y at P= 4. 

This value of y is required for the tension to increase from 
zero to half its maximum, with a corresponding change of [Gaia of 
.05 umole, so that a threshold at (Gall = .12 umole in unit of y 


equals: 


lle 


Ythreshold. .05 ‘#P ~ . 554 (a7) 


This amount can be subtracted from the equation for the activation 
factor, given by equations (2.12) or (2.13) if a threshold phenomenon 


in contraction is to be modelled. 
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2.3 Sertes Elastte Element. 

The presence of series elastic element in a muscle deter- 
mines to some extent the time course of tension change manifested 
externally. Jewell and Wilkie (1958) determined the load-extension 
characteristic of frog sartorius muscle at 0°C. Their non-linear 
load-extension curve could be divided into two regions, with va!ues 
of stiffness ranging from about 0.3 to 0.6. These are normalized 
values with respect to the maximum isometric tension and the maximum 
extension per half-sarcomere in unit of h (Section 2.1). Since the 
effect of changes of parameters in the model is studied, a standard 
value of 0.5 is given to the stiffness of the series elastic element. 
This value gave a twitch nearly identical to that computed by Julian 


(1969). 


2.4 Activatton-Contractton Coupling and Internal Shortentng. 
Activation-contraction coupling is an essential part in 
muscle function, but does not enter explicitly into the calculations. 
It is assumed that the activation process begins from the binding of 
calcium by the myofilaments, and can be separated from any previous 
events following a stimulus. Besides the threshold phenomenon, 
another difference between fast (twitch) and slow fibers in frogs is 
the presence of all-or-none action potentials in the former. The 
assumptions in the model are therefore less true for slow fibers. 
An explicit mode! for activation and activation-contraction coupling 
for these fibers might have to be included in future work. 


Although the contractile unit is kept isometric, internal 
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shortening is still possible because of the series elastic element. 
The way in which this affects the activation is an inherent property 


of the model itself. 
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CHAPTER 3 


METHODS FOR COMPUTATIONS. 

The computation method for the standard isometric twitch is 
described in this chapter. Except for the modifications given in 
Chapter 2, the method is basically similar for tetanus (repetitive 
stimulation) or when a threshold is included. A Digital Equipment 
Corporation LAB-8 computer was programmed to do the calculations. 

In the mathematical formulation of the sliding filament 
model (Huxley, 1957), the rate constants f and g for positive u are 
proportional to f; and g1, the values for f and g at u=1. For 
negative u, g = gg and is independent of u (Figure 2.2). However, 
only values for the ratio of rate constants were given in that model. 
To compute the time course of tension change, absolute values for the 
rate constants are required. From Julian's paper (1969), these are: 
f, = 81.25/sec, g,; = 18.75/sec and go = 391.9/sec. These values 
were chosen to satisfy the ratios of the rate constants and to give a 
reasonable value of the maximum unloaded shortening velocity given by 
the equation Vy, = 4mh (f; + g,), where m is the number of sarcomeres 
(Huxley, 1957, p. 288). Using an approximate value m = 10+, and 
values already given for h, fy and g, yields V,, = 4 cm/sec. 

The generator's force and equivalent stiffness given by 
equations (2.5) and (2.6) respectively can be calculated when n is 
known as a function of u. Using a small time increment At to approx- 
imate a change in the continuous time function for y, f as a function 


of both y and u, and g as a function of u, equation (2.2) becomes: 
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Equation (3.1) can be written to apply in three distinct regions over 


which u varies (see Figure 2.2): 


(i) a-e <oue<eOsd fa=s05 het] 39129 


nN = no/exp(391.9¢) 


(Ayes 0terid< VWIF Sayed (819250), 90916575 


ViCIMShoZ3) 


TES grey Pare. 7S 


VEGI IMZo) p= iV CLINI. 20) Bao. |eig 


~ [y(@) (81.25) + 18.75] exp{uly(t) (81.25) + 18.75]¢} 


Cone <n <a ern ye 10 FR a= 18.7 5u 
n(u) = no/exp[u(18.75)t] 


where t is the total time increment. 
Equation (2.5) can be written as: 
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where Pa = generator's force; A and B are as yet undefined limits. 
A simple numerical integration method is used to compute the above 
integrals. Initially a value for y is calculated at time = sAt. 
Any subsequent increment is equal to a full At. In this way the 
intermediate value is used for any interval and is assumed to 
approximate the average for that interval. The total length change, 
L, of the force generator and that of the series elastic element are equal 
because they are connected in series. This length is zero at t = 0. 
The computer then checks the particular region of u over which n 
varies and then proceeds to sum up all m in that region. In the 
regions of positive and negative u, nm is computed until a value of 
n = 0 is encountered. This occurs when nm < .0005 because of the 
word length (12 bits) for fixed point numbers in the range +1 in 
this computer. The limits-A and B in equation (3.5) are then set 
by the extent to which n is shifted into regions of negative and 
positive uw respectively. The generator's force can then be computed 
by numerical integration. The samme procedure is used to compute 
the generator's equivalent stiffness. 

At the end of a time step At, the generator exerts the 
force Pg and has an equivalent stiffness Kg. Equilibrium is established 
when the net force in the generator is equal to the force in the 


series elastic element, and is given by the following equation: 


Poe Kgl = kL OSDir) 


where k is the stiffness of the series elastic element and Z is the 
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total length change from zero time. Rearranging equation (3.7) 


gives: 


P 
eI 
ee (3.8) 
kg + k 
The force measured externally is then equal to: 
P = kL (ao) 


Addition of the stiffness of the force generator and that of the series 
elastic element results from the development by the generator in the 
active state of a stiffness which impedes both shortening and 
lengthening of itself. 

Wheel STCOS, (fetes Ci, Pg and y are stored in the 
computer. For calculating the force and length change at a later 
time, the amount of shortening or lengthening is considered when 
calculating the region for u, and the effect of m in all previous 
t is included. 

In equations (Z.12)vand C5e1) 70° (5.59, VvaluesPoteaGT, = 
5 msec and du = .025 are used. A smaller At would give a more 
accurate result, but the amount of computer time needed becomes too 
large. Halving these values does not change results appreciably 


(< 5%). The required data are printed out using a teletype machine. 
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CHAPTER 4 


RESULTS OF COMPUTATIONS. 


4.1 Abbreviattons and Definitions. 


The following abbreviations are used to describe the 


characteristics of an isometric twitch (Close, 1972}: 


tension at the peak of an isometric twitch; 
maximum rate of tension development; 
contraction time, i.e., the time from onset 
to peak of an isometric twitch; 
half-relaxation time, i.e., the time for 
tension to fall from Py to 2Pz during the 
relaxation phase of an isometric twitch; 
the time constant for the decay of tension 
during the exponential phase of relaxation 


at the end of an isometric twitch. 


4.2 Isometric Twttches. 


(t) Parameter Changes of The Acttvatton Factor. 


Figure 4.1 shows isometric twitches computed for different 
B, the rate of release of calcium from the myofilaments. 
twitch (labelled "g" in Figure 4.1) is computed using the standard 
values for the parameters given in Chapters 2 and 3, and activation 


is calculated to reach a maximum 40 msec following a single stimulus 


28 


The standard 


an RivaRiipoon cnn wide 
senyor “ita at laee it na! 


igh Met) Ors henGel nen x Heng att Fe: Mbetat 
naan ea a ore we 

hens Won) ant wt WT want? ‘nay oahas 

ee A 

mal mi? Fa | janie dbiteaatenaitad 

BI Gen yglee gaa UP gh gn! ham at wal spitht 

dite Tt yer 0F ty gona ro ghar 


tvs) een’ Season | A natin -Clieaih : 


ol eendler. 46 mabe) Let tngncipie wiciee wil pyirme 


eel bi) athe! " €o aeeane te! 


: ie 
oe 


r 


Figure 4,1 


4 8 1.2 1.6 2.0 
Time (sec) 


Computed isometric twitches for different B. The 
multiplicative factor is listed on the right of each 
curve. The curve labelled "8" is the 'standard' twitch 
for O0°C computed using 'standard' values for the 
activation factor. A tension of 1.0 indicates maximum 
isometric tension. 
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at a temperature of 0°C. It can be seen that changes in B have the 
greatest effect on the relaxation phase of the twitch, although Pt 

is also affected considerably. Figure 4.2 shows the semi-logarithmic 
plots of the fall of tension during relaxation of the different 
twitches shown in Figure 4.1. In all cases the tension record shows 
a final exponential fall. Jewell and Wilkie (1960) observed this 
exponential fall of tension in relaxing muscle, and gave evidence 
that the rate constant of tension fal! appeared to be characteristic 
of a particular muscle under given conditions. The rate constant (1/7) 
of fall of tension from 60% to 5% of Py of the standard twitch is 
about 5.75/sec. For twitches computed with different 8, these rates 
are normalized to the standard rate and are listed in Figure 4.2. 

The relaxation rate is approximately doubled with doubling of 8 
except for the highestrate (8 x 4). 

The effects of other activation parameters have also been 
computed. These include Co, the amount of calcium released by one 
stimulus and a, the rate of binding of calcium to the myofilaments. 

A fourth one, designated SR, stands for volume of sarcoplasmic 
reticulum, assumes that both Cg and 8 are changed by the same factor. 
It is assumed here that the existence of a smaller volume of sarco- 
plasmic reticulum will cause a smaller amount of calcium to be 
released, from both the sarcoplasmic reticulum itself and the bound 
state of the myofilaments. In the computation, both Cg and 8 are 
changed arbitrariiy by the same factor. 

To compare the effects of parameters changes of the acti- 


vation factor on the characteristics of isometric twitch, the para- 
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Figure 4,2 


0.1 


Tension (P) 


0: ne) 20 30 
Time (sec) 


Semi-logarithmic plots of the fall of tension during 


relaxation in the computed twitches shown in Figure 4.1. 


The relative rate constants of fall of tension from 60% 
to 5% of the respective maxima are: 3.2 (B x 4), 2.2 
CBZ) sl Boi. aol ap eee pees) Ope aca) « 
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meters and the characteristics relative to their 'standard' values 


are plotted on a log-log scale as shown in Figure 4.3. Parameters 


of the activation factor changed by a factor of two and four are 


shown. 


(it) 


(1) 


C2) 


(3) 


(4) 


Several results follow immediately from Figure 4.3: 


The largest change observed is the change in Pin 
for B < 1, where a direct proportionality is ob- 
served. A proportionality was also observed for 

the rate constants of tension fall with different 

8, as previously mentioned. The other characteristics 
appear to depend on several parameters to a varying 
extent. 

The time course of the contraction phase as defined 

by dP/dt and Te is more sensitive to changes in Cg 
OGeoh. 

The time course of relaxation as defined by Tin 
is most sensitive to change in 8. 


a seems to have the least effect on the characteristics 


of the tsometric twitch. 


Parameter Changes of the Original Sliding Filament Model. 


Figure 4.4 shows results obtained by changing independently 


the parameters in the original sliding filament mode! (Huxley, 1957). 


These are the boundary rate constants for the making and breaking of 


cross~bridges. 


Also included is the stiffness of the series elastic 


element which does not exist in the original model. Its effect wil! 


be considered 


in more detail in the next section. Several obvious 
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Peak twitch tension (P,) 
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Max. rate of tension development(dP/dt) 


Figure 4,3 
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Relative value of parameter 


Effects of changing parameters in the activation factor 
on the characteristics of the isometric twitch. All| 
values have been normalized to their respective 'standard! 
values. 
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Figure 4.4 
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Relative value of parameter 


Effects of changing parameters in the original sliding 
filament model and k (series stiffness) on the char- 
acteristics of the isometric twitch. AIl values have 
been normalized to the 'standard! values. 


5 — ah 


7 


- 7 
a 


i 


* aa 


gate qeeneaenns FEAT) 


results are outlined below: 

(1) Similar to the parameters of the activation factor, 
the twitch characteristics depend on several para- 
meters and are affected less than proportionally by 
changes in any one of them; 

(2) The peak twitch tension is rather insensitive to 
any one parameter; 

(3) The contraction phase is more sensitive to fy or k. 

(4) The relaxation phase is most sensitive to g}. 

(5) go has negligible effect on the characteristics of 


an isometric twitch. 


(tit) Effect of the Sertes Elastic Element. 

Figure 4.5 shows the isometric twitches computed for 
different kK. In the model k has the greatest effect on the con- 
traction phase (see also Figure 4.4). Jewell and Wilkie (1960) 
reported that k had very little effect on the exponential relaxation 
phase. Hill (1951a) showed that an added series compliance to an 
‘isometric! twitch of frog sartorius affected the time course of 
tension development. To compare, the series stiffness k in the 
model was converted to absolute units of compliance. Added compliance 
was obtained after subtracting the compliance of the standard twitch 
which was assumed to have zero added compliance. Its effect on the 
characteristics Py, dP/dt and T, is shown in Figure 4.6. Data from 
Hill1951a, Figure 1) on the effect of an added compliance on the 


twitch are also plotted in Figure 4.6. The standard value of .5 for 
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Figure 4.5 Computed isometric twitches for different k. The 
multiplicative factor is listed with each curve. The 
curve labelled "k" is the 'standard' twitch. 
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Figure 4.6 
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Effects of added compliance on the characteristics of 
the isometric twitch. 0, experimental data from Hil 
(1951a); @, data from model. 
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k is equivalent to a compliance of .018 mm/gm. Good agreement is 
found between the two sets of data for P+ and dP/dt. However, the 
series elastic element in the model produces greater changes in 7, 


than are found experimentally (Figure 4.6C). 


4.3 Tetanus (Repetitive Stimulation) and Threshold. 

Figure 4.7 shows the tension curves obtained by stimulating 
the 'standard' muscle at a rate of 10/sec, together with a muscle 
with zSR as defined previously. I+ can be seen that the tetanic tension 
curves follow the same time course initially as the corresponding 
twitch up to .1 sec. At the end of each stimulus, Ap and a new 
value of Co,w are computed. These new values are then used to com- 
pute the tension curve for the next stimulus. The fraction of 
interacting cross-bridges due to the previous stimulus is retained 
by the computer. As explained in Chapter 2, Co, Ms the amount of 
calcium released by one stimulus, becomes smaller with each additicnal 
stimulus. 

eer differences are obvious between the two tension 
curves; the muscle with 4SR attains a smaller twitch tension, has a 
longer contraction time and a lower rate of tension development, but 
the tension during repetitive stimulation is better maintained. 

As defined by Julian (1969), the largest value y could 
attain in the 'standard' activation factor is 1, but with repetitive 
stimulation in the present model, y can be several times greater. 

From equation (2.7), this means that f will also be increased several 


times. However, the generator's force can only reach a certain limit 
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Figure 4.7 
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Computed isometric twitch and tetanus (10/sec) of 
A) 'standard' muscle and B) 'sSR' muscle. The 
peaks of the tetani are also marked. 
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for the following reason. From equation (2.2), it can be seen that 
as f becomes very large, n approaches a limiting value of 1. Sub- 
Stluiuringi2 = i and (f + ¢)/f = .8125 into equation (2.5) end eval— 
uating Pg using the limits of integration for isometric steady-state 


gives: 
1 
Pg = 2(.8125) fo udu = 1.23 (4.1) 


This is the maximum force which the generator can attain when there is 
no internal shortening. The force P measured externally is, however, 
always smaller than Pgs as is obvious from the following equation 


after substituting L from equation (3.8) into equation (3.9): 
ie (4.2) 


It is also obvious that when k is very large (approaching a true iso- 
metric condition) P > Po: 

Indirect experimental evidence on human muscle indicates 
that the active state cannot only be prolonged, but can also be 
intensified (Desmedt and Hainaut, 1968). 

The activation factor is calculated for 0°C. At this low 
temperature the twitch/tetanus ratio is high and equal to .81 for the 
"standard! muscle compared to a !ower value of .70 for the 'zSR'. 

A threshold value of .554 for y is derived in Chapter 2. 

If this is deducted from the activation factor for 'sSR', no tension 


will be generated when only a single stimulus is delivered. This is 


wee ee " me 2 Rh pi paite re ers 
CPs at Jr “6 a = ifort ae 
- 


Lae 
oe 
an me! ina fetieea yal ' a 
Gutre~ 4 hi ers cde ~a fe jrencatnl eo 


— 
7 uw ' , 
‘sre sii arc te mm Qf 1 restate aft et oe 


@ | oi% : ry S - 
Part) 6 « 7 ae ner on 
q | we ge | 4 oi eR ’ et : , wri J +4 ° - ’ ral : 
, 
ar ers ry 7 an wet la) @ i ae a1 ¢ 2) =“ ot Ter! el (oem zy ol 
J a : 2.78) iP id ape nee 
- ~ 2 - a _ = — 
. . i) a.) ee 
a 9 is 
y , 
, 5‘ 4? ——- Sa 8 = 4 
T » > 7» - > 
i, - 7 
i ‘. 
n)} ‘ ; - 
—fe) MAN OTP itoeovagal, errr yes St ree owt ‘gale Veo oete e o 
é glee ; 
: —» = seit th i, 20 
- . re a , hs ii it 
walhe!| eto4u) hanes? _ > S5rta! /2 leteees vay * Bas Hew i : 


66 Api neo Tye yhepncie: 20,8 al ywine t¢ med ours sd i bit 


oa”: ph) Pa e-* ad) = = 


daw tweets ane > eee ai) A 


qi 
coe 
_ “- 


* = . ie 


ay 7 pat’ FD jor ois Ohad) oi wep hey ai sah th ca 
¥ a “ 
ue i aliens bak ay ryt 


aft a ita ree : 
Gs my ,’ er Salles tna i | 

mane es oer Pt. ko thd 
eae ine stem we ou _— 


i x6 ae ay ws 
ot eld? era 1 a) ay 


4] 


because Cg,; is reduced by half in the calculation of y in'4SR' and 
- the maximum value for y in this case is .5, which is smaller than 
.954. When this is encountered, y is set equal to zero in the com- 
putation. However, in repetitive stimulation y can exceed the 
threshold value of .554 and tension can then be generated. Figure 
4.8 shows the tension curve of a 'standard' and a '$SR plus 
threshold' muscle stimulated at a rate of 20/sec. The results are 
similar to that at 10/sec although at higher rates the effects are 
more prominent. 

The muscle with 'sSR plus threshold' does not generate any 
tension until after the second stimulus. Comparing with the 'standard' 
muscle, it also reaches a smaller peak tension at a longer time, a 
smaller rate of tension development, a smaller fusion frequency and a 


low fatiguability. 


4.4 Frequency-Tenston Characteristic. 

Muscle tension increased with increasing frequency of 
stimulation in a characteristic way (Cooper and Eccles, 1930; 
Matthews, 1959; Rack and Westbury, 1969). The frequency-tension 
curve has a sigmoid form. There is always a middle range within 
which tension rises steeply with each increment in frequency, but 
further increase in stimulation beyond that range has diminishing 
effects on the tension. 

lf the parametric modelling of activation for the three 
types, 'standard', 'zSR', and 'zSR plus threshold', represents 


different patterns of activation ina muscle, or in different muscles, 
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Figure 4.8 
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Computed tetanus (20/sec) of the 'standard' muscle and 
the 'sSR + threshold' muscle. The peaks of the tetani 
are also marked. Note tetanic tension in the '4S5R + 
threshold! muscle begins to develop only after the 
second stimulus. 
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their frequency-tension relations will reveal the characteristic 
pertinent to each type. At all frequencies of stimulation the 
tension usually reaches a peak and then declines, whether the 
tetanus is partially fused or completely fused. This is due to 
the residual Ag and n from the previous stimulus and a proportion- 
ally smaller Cg for each additional stimulus. For this reason, 
the portion of tension record containing the peak and between two 
stimuli is averaged with respect to time for the partial ly-fused 
or completely-fused tetanus, while in a twitch the time-average of 
the whole record is considered. 

Figure 4.9 shows the frequency-tension relations for these 
three types of parametric modelling obtained in this way. The curve 
for the 'sSR plus threshold! has a sigmoid shape while the curves for 
the 'standard' and '4S5R' rise very sharply. 

In real muscles the sigmoid curve is always observed for 
the frequency-tension relationship. This is true in a slow twitch 
muscle such as the soleus muscle (Matthews, 1959; Rack and Westbury, 
1969; Mannard and Stein, 1972), or in fast twitch muscle such as 
the gastrocnemius in cat (Kernell, 1966). The maximum rate of rise 
of the frequency-tension curve is greater in a fast muscle and 
Matthews (1959) reported by stretching the soleus muscle, the curve 
shifted towards the lower frequency end (the muscle had a lower fusion 
frequency) and had a greater maximum slope. In the present mode! the 
existence of the sigmoid shape is a result of the threshold. 

Tension does not appear until frequency is just above 1!/sec 


in the 'ISR plus threshold' muscle for the reason which is given earlier. 
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> SR+ Threshold 
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Frequency of stimulation( cycle per sec) 


Frequency-tension curves of three muscles with different 
activation mechanisms: @, 'standard'; 6, '4SR'; A, '45R 
+ threshold’. Note sigmoid shape of curve for the 

'lSR + threshold’ muscle and absence of tension below 


1 cycle/sec. 


4 


oe Set 


a! ae re gil 
ae ek wh 
4 


| ye hi | 
7 “ bee 


‘a ee 
\ ot 
Nyt heel A 


45 


At low frequencies, the 'SR' muscle actually gives a 
greater average tension than the 'standard' muscle. This is because 
although Co in the '3SR' muscle is only half of that in the 'standard' 
muscle, the smaller 8B greatly prolongs the duration of activation and 
thus the time course of the relaxation phase. Therefore, summation 
of force occurs from one stimulus to the next. At high frequencies 
the two curves are almost identical; here the peak tension is 
approaching the limiting value (1.23) of the maximum possible tension 
in the model. 

In the soleus muscle of the cat at 37°C, a stimulus frequency 
of 30/sec already gives a tension more than 90% of the maximum at the 
optimum length (Matthews, 1959; Rack and Westbury, 1969). In the mode! 
a seeposnerd ihe frequency would be about 10/sec. This is mainly 
because activation is originally derived for 0°C, a temperature at 
which fusion frequency is low and activation reaches its peak even 
in a single twitch, and partly because the greatly increased y tends 
to make n approach its maximum value of 1, thus giving a high tension 


even at a frequency of 10/sec. 


4.5 Fourter Transform of the Isometric Twttch. 

The direct Fourier transform is a tool that resolves a 
signal into its frequency components, and gives a frequency domain 
representation of the signal. Time domain representation specifies 
a function at each instant of time, whereas frequency domain repres- 


entation specifies the relative amplitudes of the frequency components 
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In a nerve-muscle preparation, the action potentials can be 
considered as a series of unit impulses (French and Holden, 1971a; 
Knox, 1969) and as an input to the muscle which develops tension as 
its output. Since a motoneuron usually conveys information in the 
form of its impulse frequencies, frequency domain analysis applied 
to the neuromuscular system is of great value in revealing the signa!- 
handling characteristic of a muscle (Partridge, 1966). Methods have 
been developed by French and Holden (1971la, b, c) for estimating the 
frequency response of a muscle using spectral analysis. Alternatively, 
if the nerve-muscle preparation were considered as a linear system, 
the twitch could be used as the impulse response of the system, and 
the Fourier transform of the twitch would give the frequency response 
(D'Azzo and Houpis, 1966). 

A Fourier transform generates complex numbers which need two 
plots for their graphical representation as a. function of frequency, e.g. 
plots of magnitude and of phase. Figure 4.10 shows the magnitude of The 
Fourter transformand the phase angle as a function of log frequency 
of the standard isometric twitch. The magnitude plot essentially 
gives a prediction for the responsiveness of the muscle in terms of 
the change in tension which should be produced for each impulse per 
second modulation in stimulus rate, if the latter were being modulated 
sinusoidally with the frequency given on the abscissa. The fitted 
curves in Figure 4.10 are the magnitude and phase which would be 


expected for a critically-damped, second-order system of the form: 


G(s) = (4.3) 
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Figure 4.10 The magnitude (A) and phase (B) of the frequency response 
of the 'standard' muscle obtained by Fourier transformation 
of the isometric twitch. Further explanation in text. 
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where s is the Laplace transform variable; G(s), the transfer function 
representing the frequency response; Go, the zero frequency gain; f,, 
the undamped natural frequency; and z, the damping ratio of the system. 
GoGechitical damping, ~iis equal to.1. 

A second-order system is completely defined by the three 
parameters Go, f, and t in its gain curve which is a plot of the mag- 
nitude of G(s) against frequency (the continuous curve in Figure 4.10A). 

The zero frequency gain is the magnitude of the direct 
Fourier transform when frequency approaches zero and determines the 
vertical position of the curve. It is equal to the area under the 
tension record of a twitch as can be seen from the following equation: 


~2nft 


lim F(f) = lim S72 F(t) e dt = 5°" F(t)dt (4.4) 


fo fro 


where f(t) is the tension curve (in the time domain) and F(f) is 
the Fourier transform. 

The undamped natural frequency is the frequency with which 
the system will oscillate if ~t = 0 and determines the horizontal 
position of the gain curve. It Is the frequency at which Go 
decreases to half its value in the gain curve for a critically- 
damped, second-order system. It is related to the speed of response 
of the system in the time domain. Two characteristics which deter- 
mine the speed Sper aven change in an isometric twitch are the max- 
imum rate of tension development and the rate of tension decay. The 
larger the fy,, the faster the system responds. 


The damping ratio while specifying the shape of the gain 
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curve in the frequency domain, also determines the shape of the twitch 
in the time domain. A critically-damped system has t = 1; any value 
smaller or greater than 1 indicates, respectively, underdamping or 
overdamping. In a second-order system given by equation (4.3), the 
impulse response (the inverse Fourier transform) for ¢ > 1 is highly 
skewed whereas for t < 1 some overshoot (oscillation) is observed 
Ogata, 1970); 

If the frequency response of the standard muscle is rep- 
resented as a critically-damped, second-order system, it has a 
natural frequency of about 1 Hz (Figure 4.10A). The steep fa! l-off 
of the Fourier transformat high frequencies effectively describes 
also the low-pass characteristic of the contractile unit in the model. 
In réal muscles this characteristic has been reported in frequency 
response studies under isometric conditions (Partridge, 1965; 
Poppele and Terzuolo, 1968; Rosenthal, McKean, Roberts and Terzuolo, 
1970; Mannard and Stein, 1972) and isotonic conditions (Partridge, 
1966). 

Figure 4.11 shows the Fourier transforms of isometric 
twitches computed by changing 8, SR and k by a factor of two. The 
response of the standard muscle is included for comparison. The 
representation by a critically-damped, second-order system is also 
good in most of the Fourier transforms shown, although in some cases 
a curve for a slightly underdamped (ct = .9) or a slightly overdamped 
(c = 1.2) second-order system gives a better fit. The second-order 
system parameters of the magnitude plots shown in Figure 4.11 are 


listed tn Table 4.1. 
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Figure 4.11 The magnitudes of the Fourier transforms of computed 
isometric twitches after changing different parameters: 
A) 8, B) Sk and C) ko +, 'standard!: A, doubliing of 
parameters; 0, halving of parameter. 
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TABLE 4.1 Comparison of the second-order parameters used to describe 
the frequency response obtained by Fourier transformation 
of the computed isometric twitches after changing differ- 


ent parameters by a factor of two. 


Parameter Relative Go fn (Hz) Estimated f 
(at .1 Hz) 

Standard 0 120 Tsu 

(Spe one OReke 
B x 2 1.88 0.56 ape 
Boxe? O254 en) 0.9 

Sk x + Pee 0.60 122 

Sh xX Z O75 1.50 0.9 
kx 0.94 0.95 Oe) 


ays 


In the twitch tension records for different 8 (Figure 4.1) 
and from Figure 4.3, it can also be observed that the characteristics 
of the twitch are most sensitive to changes in B or SR, which change 
the characteristics in nearly the same fashion. Decreasing 8 or SR 
produces increase in the 'sluggishness' (increases in teanc Tin 
in general) of the muscle. 

The Fourier transforms are also more sensitive to 8 and SR 
than to k and are similar in that both 8 and fy, change in the same 
fashion. 6 has a greater effect because its 'standard' value is 
smal! (6/sec), and changing it by a factor of two has a relatively 
greater effect on the isometric twitch, and hence the Fourier trans- 
form. In the case of SR,Cg is also changed by the same amount, and 
this tends to reduce the maximum value of y and hence the effect of 
B on the time course of y. The increased 'sluggishness' as a result 
of decreasing 8 or SR produces concomitant increase in t in general 
and decrease in f,. Gg, which measures the area under a twitch 
record, depends not only on Tg and Tip, but also on Py. 

As shown in Figure 4.11C, changing k, the stiffness of the 
series elastic element, by a factor of two has only a minor effect 
on the Fourier transform. An increase in k has negligible effect on 
t, although Gg and fy, are slightly increased. Recent studies on the 
frequency response of isometric soleus muscle in the cat (Mannard and 
Stein, 1972), however, showed that increasing the stiffness by 
stretching the muscle increased Go and ct, while f, was decreased. 

As shown in Figure 4.5, increasing k results in a greater dP/dt and 


smaller T,; these changes are consistent with an increased fy, 
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observed for the gain curve. Changes in the form of the twitch and the 
area under the tension curve are probably too small to have a signifi- 


cant effect on the Fourier transform. 
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CHAPTER 5 
METHODS FOR EXPERIMENTATION. 


5.1 Preparatton. 

All experiments were performed in sttu on the sartorius 
muscle of the bullfrog, Rana catesbiana. Preliminary experiments 
showed that if circulation was kept intact, rapid fatigue could be 
prevented. 

The animal was anesthetized by immersion in a .2% solution 
of Tricaine Methanesulfonate (Fraser, Vancouver) for 30 to 40 minutes. 
Narcosis lasted several hours. The distal tendon of the right sar- 
torius muscle was exposed and freed of the adjacent connective tissues. 
An inextensible steel hochewas tied directly to the distal tendon. 
Particular attention was given to tying the silk thread as tightly 
as possible very near the ends of the muscle fibers, thereby reducing 
the series compliance. 

The right sciatic nerve distal to the pelvis was exposed 
from the dorsal surface. The large collateral branch which ultimately 
innervated the sartorius muscle was cut together with a few millimeters 
of sciatic nerve. All other muscles of the right limb were denervated. 
Care was taken to avoid puncturing any major blood vessels. The 
distal tendon was cut with the steel hook attached to it. 

The animal was put into a Perspex trough in its prone 
position. The trunk was secured by pins forced against its pelvic 


bone and the legs were clamped to one side of the trough. In this 
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position, the right thigh was at an angle of approximately 150° from 
the side of the body and the knee was about 1 cm above the floor of 
the trough. The thigh muscles near the dorsal surface were spread 
apart to make a paraffin oil bath in which the nerve was stimulated. 
Ringer's solution (NaCl, 117 mM; CaClo, 1.0 mM; KjHPO,, 2.4 mM; KH»PO,, 
0.6 mM; 200 mg% glucose) containing .03% of the anesthetic Tricaine 
Methanesulfonate was poured into the trough so that the animal was 
about half immersed. Under this condition the animal could be kept 
unconscious indefinitely. Limb movements were absent but blood cir- 
culation in small cutaneous vessels could be observed. A mixture of 
95% Op and 5% COs was continuously bubbled into the Ringer's solution. Temper- 
ature changes between 3° and 31°C were achieved by circulating water 
in glass tubings submerged in the bath. The temperature of 

the circulating water was regulated by a thermistor ee TOMa 
cooler (Precision Scientific Lo-Temptro! 154). Temperature of the 
Ringer was measured to the nearest half degree. 

A muscle was stretched from a slack length until it gave a 
maximum twitch tension. This length (Lo) was not changed throughout 
an experiment. The weight and Zo of the eight muscles used in the 
experiments were 0.97 + 0.07 gm (mean + S.D.) and 6.55 + 0.44 cm, 


respectively. 


SO LLIMULOCLON. 


Two types of electrical stimulation of the muscle nerve were 
used in the experiments. Regular pulse trains were generated by two 


sets of digital counters and stimulators (Bell and Stein, 1971). 
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Stimulation using one frequency could be switched to another almost 
instantly. A random pulse train was produced by adding random noise 
to a neural analog (Stein and French, 1970) which produced regular 
pulses in the absence of the random input. These pulses were used to 
trigger a stimulator. The random pulse train used had a mean rate of 
5/sec and the intervals between pulses range from about 25 msec to 

1 sec. Stimulus pulses of twice maximal strength and 100 usec 
duration were used. The stimulating electrodes were made of a pair 


of pure silver wires 0.5 mm in diameter. 


5.3 Recording and Analysts. 

lsometric tension was recorded using a Grass FT-03C trans- 
ducer and a Tektronix 3C66 carrier amplifier. The transducer had a 
compliance of 1 u/gm and a resonant frequency of about 330 Hz. Tension 
records were simultaneously examined on a Tektronix RM565 Dual Beam 
oscilloscope, monitored on a Hewlett-Packard 7700 pen recorder and 
stored on magnetic tape (Thermionic Products, T-3000 FM recorder) for 
analysis. 

A Nihon Kohden PC-2A camera was used to photograph isometric 
twitches in sequence from the oscilloscope face. Peak tension, con- 
traction time and half-relaxation time were measured by hand. 
Frequency domain analysis was carried out using a LAB-8 computer 
(Digital Equipment Corporation). Spectral analysis packages (French 
and Holden, 1971a, b, c) were available to evaluate the frequency 
response given any analog or digital input and output. In the 


analysis here the input and output are respectively the random pulse 
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train and the tension generated. 


5.4 Descriptions of the Frequency Response. 

The parameters used to specify the gain curve of the frequency 
response of a second-order system are given in Chapter 4. These are 
zero frequency gain (Go), natural frequency (f;,) and damping ratio (z). 

The data points for the estimated frequency response were 
converted to units (gm-sec)/impulse at a certain frequency, displayed 
on the screen of a storage unit (Tektronix 601) and photographed using 
a Polaroid camera. 

A family of 16 frequency response curves for the gain of a 
linear, second-order system was constructed on a transparency with 
logarithmic axes having the same scales as those of the Polaroid photo- 
graphs. By sliding the transparency over the photographs, the curve 
providing the best fit to the gain data was selected and the parameters 
Go» fy and ¢ could be specified to an accuracy of about Pepa 

A coherence function which gives a normalized measure of the 
extent to which the linear frequency response function accounts for the 
total response of the system at any given frequency was also printed out. 
A coherence of 1 would indicate the behaviour of the muscle could be 
described by a linear, time-invariant system.with no noise or other 
extraneous inputs. The data were also printed out using a teletype 


machine. 
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CHAPTER 6 


RESULTS OF EXPERIMENTATION. 


6.1 Effects of Tetanus and Temperature on the Characteristics of the 
Isometrte Twitch. 
(t) Peak Twttceh Tenston. 

The effect of tetanus on the isometric twitch tension obtained 
by stimulating the muscle nerve was studied in the temperature range 
3z to 30°C. A tetanus of a few seconds duration produced a potentiation 
of the isometric twitch (a greater Pz and a longer time course), but the 
amount of potentiation was not maximal until some time after the end of 
the tetanus. Figure 6.1 shows the isometric twitch tension record of 
one typical experiment at 20°C. Immediately after the tetanus, peak 
tension rose nearly linearly to a peak and decayed with an approximately 
exponential time course. This potentiation after a delay was also 
observed in frog sartorius by Connolly, Gough and Winegrad (1971), when 
the muscle was stimulated directly. In the record shown in Figure 6.1, 
the rising phase of potentiation lasted about 16 seconds, and was 
followed by a falling phase that lasted for over 3 minutes. 

The effect of tetanus produced by about 250 stimuli was 
observed at any one temperature. At temperatures above 10°C, a 
stimulation rate of 50/sec was used, while at lower temperatures, 
33/sec was used. These stimulation rates were adequate and satisfactory 
to produce potentiation in these temperature ranges. Close and Hoh 


(1968b) showed that the total number of stimuli ts the more important 
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Figure 6.1 The effect of tetanus on the isometric twitch tension 
of a sartorius muscle at 20°C. Rate of stimulation was 
30/min following a 5 sec tetanus at 50/sec. 
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parameter in determining the degree of potentiation a rat fast 
skeletal muscle. Preliminary studies indicated that this was also 
true in frog muscle fibers, although at lower temperatures (below 
10°C), a lower rate of stimulation (33/sec) seemed to be more 
effective. 

To produce twitches, the nerve was stimulated at a rate 
of 30/min. At lower temperatures, even this rate would somewhat 
alter the characteristics of the twitch; therefore, only about 10 
twitches were elicited before the tetanus. The isometric twitch 
immediately before the tetanus was referred to as the ‘control! 
twitch. After a temperature change, about 30 min were allowed for 
the muscle to equilibrate with the new temperature. 

The separation of the post-tetanic potentiation into two 
phases was a common feature for a wide range of temperatures, although 
at temperatures greater than 25°C the rising phase was usually absent. 
This means that maximal potentiation was observed in the first few 
seconds after a tetanus. In the case where a rising phase was 
present, the value for peak tension at 2 sec after a tetanus was not 
always the same as the pretetanic value; sometimes it was higher and 
occasionally lower. 

Figure 6.2A shows, in a linear plot, peak tension of the 
isometric twitches in the rising phase at time t after the tetanus 
at different temperatures for one muscle. The points at any one 
temperature are fitted by straight lines. P, of the Leonaiods 
twitches are also shown. The peak tensions of the maximallypotentiated 


twitches are shown separately in Figure 6.4. 
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Figure 6.2 A. Linear plots of time course of rising phase of poten- 
tiation at different temperatures. Peak tension, Pz, is 
shown at time t after a tetanus. 'Control' indicates pre- 
tetanic values. The last point of each set has a value 
about 80% of the maximum value attained. 

B. Semi-logarithmic plots of time course of falling 

phase of potentiation. Difference in peak tension is 

shown at time t after a tetanus. The first and last points 
of each set have respective values of about 80% and 30% 

of the maximum difference in Py. Temperatures: 4° (A), 

OF°8 (0), lg? (0), 24> (yo vand 292-6 (7) Oniy Iwepot iis 
are shown for the falling phase at 4°C because of the min- 
imum difference (1.25 g) distinguishable in the measurement; 
two other points were also used to obtain the fitted line. 
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Figure 6.2B shows, in a semi-log plot, the difference in 
Pz of the twitches during the falling phase. The difference in Py 
was usually obtained by subtracting Pz of the 'control' twitch from 
that of the twitches in the falling phase. At higher temperatures, 
it was not unusual for the post-tetanic Pz to fall to a lower value 
after the potentiation and increase again, or remain at that level. 
In this case, the difference was obtained by subtracting the Peein 
the 'depressed' twitch. The differences of peak tension obtained in 
this way for the twitches in the falling phase are shown in Figure 
6.2B. The time scale indicates the time course of the fall ¢ sec 
after the tetanus. 

Different rate constants (in units of 'gm/sec' for the 
rising phase and 'sec !' for the falling phase) were obtained from 
the fitted lines on the two phases. The rate constants of change 
of the rising phase of potentiation at different temperatures for 
four muscles are shown in Figure 6.3A, in a semi-log scale. Temp- 
erature changes for the four preparations were as follows: increasing 
temperature in two muscles, decreasing in one, and increase £61 lowed 
by a decrease in the fourth. Data were obtained only for temperatures 
below 25°C because at higher temperatures, maximal potentiation was 
usually reached in a few seconds after a tetanus and too few twitches 
wereobtained to give an accurate measure of the rising phase. Below 
25°C, the rising phase of potentiation is fitted by a straight line 
(corr. coett. =~:9/)eand NaS.a GigeCneaboun 2aJelo.D eae 17.2 

Figure 6.3B shows the rate constants of differences in Pz 


of twitches in the falling phase in the same experiments. The points 
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Figure 6.3 
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A. Semi-logarithmic plot of rate of increase of peak 
tension (rising phase of potentiation) at different 
temperatures. 

B. Semi-logarithmic plot of rate constant of decay 

of difference of peak tension (falling phase of poten- 
tiation) at different temperatures. Data are obtained 
from the fitted lines shown in Figure 6.2 and from 
three other experiments. 
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are also fitted well by a straight line (corr. coeff. = .98); a 
Qyp OF 239°(5.D. = .18) [s= obtained from the fitted line. 
Figure 6.4 compares the peak tension of the 'control' and 
the maximallypotentiated twitch at different temperatures, obtained 
from the same experiments. The two sets of points are fitted by 
straight lines. Since the weight and optimum length of any of the muscles 
used deviated within 10% of their mean values (see Chapter 5, Section 
5.1), and the error of measurements in the experiments was estimated 
to be about +15%, the factors of weight and length were not considered. 
Several results are obvious from Figure 6.4. 
(1). Tension varies over a wide range of temperature. 
(2) Peak tension decreases with temperature. 
(3) For the 'control' (pretetanic) twitch, peak tension 
decreases by about 49% for a temperature change from 
SY dg oes Op Cp 

(4) For the same temperature change, peak tension of 
the maximally potentiated twitch decreases by about 
22%. 

The decrease of P+ with temperature was observed for frog 
musclessinatne range O- to 20°C CHINT, 1951cs Kelley andeiry, 1953) 
whereas Jewell and Wilkie (1958) reported both an increase and 


decrease in Pz. 


(ii) Contraction Time, Te, and Half-Relaxation Time, Ty. 
The contraction time and half-relaxation time of the 


‘control! twitch were greatly decreased with an increase in temperature. 
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Figure 6.4 
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Peak isometric tension (q ) 
m 
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The relation between temperature and peak tension in 

the 'contro!l' twitch (0) and the maximally potentiated 
twitch (@). Each set of points is fitted by straight 
lines (least square method). Tension decreases at 
Aubate ot) 14.9) 0/1050 fot ss.0mc ar sOcGa ine the sscon trom 
twitch (---) and 7.6 g/10°C from 93.7 g at O°C in the 
maximally potentiated twitch (——). Data from four 
sartorius muscles. 
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Tetanus also produced an immediate increase in both Jo and Tin. 
The values then returned very slowly with a variable time course 
to their control values. 

Since Tf, and Tin represent the duration of tension change 
of the contraction and relaxation phase respectively, their reci- 
procals indicate rates of change. The reciprocals 1/T@ and 1/TLy 
were referred to as the contraction rate and relaxation rate res- 
pectively and were plotted in Figure 6.5 at different temperatures 
for the 'control' twitch and the maximally potentiated twitch from 
the same experiments as in the previous section. 

Each set of data points is fitted by straight lines (all 
have corr. coeff. > .9). If semi-logarithmic fit were used instead, 
the correlation coefficients are not significantly different and 
the @19's so obtained areuinetne rangeal.9 to" 2.2." Using rhe slopes 
of the fitted lines, the ratio T1,/T, for the whole temperature range 
is found to be equal to 1.3 for the 'control' twitch whereas a value 
of 2.0 is obtained for the maximally potentiated twitch. 

Using the values from the fitted lines, it can be observed 
that the ratio May th is increased by tetanus from 1.8 to 1.9 at 
10 Gy and trom 1.5 to 2.2:at 25°C. An other words; although both 
Tip and Te are increased by the tetanus, the percentage increase in 


Tin/Tg is greater at higher temperatures. 


6.2 Effects of Temperature and Tetanus on the Frequency 
Response of the Isometrte Muscle. 


The response of the isometric muscle to nerve stimulation was 
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/T, (sec ) 


Temperature (°C) 


The relation between temperature and A) contraction rate, 
1/fo2 and B) relaxation rate, 1/Tiy in the ‘control! twitch 
(--0--) and the maximally potentiated twitch (—@—). 


Data from four sartorius muscles. 
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studied in the frequency domain. Figure 6.6A shows an isometric 
twitch of sartorius obtained by supramaximal stimuiation of the 
muscle nerve at 20°C. The magnitude of the Fourier transform com- 
puted from this twitch by the method of French and Holden (1971a, 
b, c) is shown in Figure 6.6B. 

The fitted curve in Figure 6.6B is the response which 
would be expected for a critically-damped, second-order system which 
has a transfer function given previously by equation (4.3) in 
Chapter 4. The critically-damped, second-order representation is 
extremely good at this temperature, although in some cases a4 curve 
for a slightly overdamped second-order system gives a better fit. 
Underdamping is not observed in the frequency response of real muscles. 
The frequency response has a shape which resembles the Fourier trans- 
form of the 'standard' twitch described in the model previously. 

Figure 6.6C compares the observed phase angles with those 
derived from the frequency response curve in Figure 6.6B. The 
deviation from the predicted angle results from the delayed response 
of the muscle to stimulation. Time delay does not affect the mag- 
nitude of the frequency response, but increases the phase angle in 
proportion to its frequencies (Ogata, 1970). However, when the ob- 
served 6.5 msec latency between muscle nerve stimulus and muscle 
response is taken into account, the dashed curve in Figure 6.6C is 
obtained. This indicates that if the stimulus is applied to a model 
with the fitted frequency response, the output in the time domain 
closely resembles the tension record obtained experimentally. The 


latency in the experiment results mainly from nerve conduction, 
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Figure 6.6 
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aon 


100 msec 


Tension 50g 


0.05 


-180 


Phase ( deqrees) 


-270 


Frequency ( Hz) 


lsometric twitch (A) of sartorius muscle obtained by 
stimulation of the muscle nerve at 20°C. The magnitude 
(B) and phase (C) of the frequency response of sartorius 
muscle obtained by Fourier transformation of the iso- 
metric twitch. Further explanation in text. 
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synaptic transmission and Reiter ones coupling. 

As pointed out earlier, the frequency response computed 
from the isometric twitch is adequate to describe the responsiveness 
of the muscle to nerve stimulation if the system is linear. In the 
nervous and muscular systems, many non-linearities exist. Moreover, 
nerve impulse trains in motor axons are variable in time and asynchronous. 
This can be approximated to some extent by applying random stimulation 
to the muscle nerve and recording the partially-fused or unfused 
tension. 

The effect of tetanus on the frequency response was deter- 
mined using this method in the temperature range 3° to 31°C. The 
stimulation sequence was as follows. A random pulse train that had 
a mean rate of 5/sec was applied to the muscle nerve for 30 sec and 
tension was recorded simultaneously (referred to as the 'contro|' 
record). At least 15 min were then allowed to elapse. This permitted 
the effect of the random stimulation to subside and the muscle to 
return to its resting condition. A second random train was then 
applied immediately after a tetanus elicited by about 250 stimuli in 
a way given in Section 6.1; the tension recorded was referred to as 
the 'post-tetanic'! record. In all the records used for analysis, 
the peak tension of the isometric twitch before the first random 
stimulation and that before the tetanus did not differ by more than 
20%. Temperature change was achieved in a manner given before. 

Figure 6.7 shows the frequency responses obtained from a 
‘control! and a 'post-tetanic’ record at 16°C. Tetanus produces an 


increase in Gp by about 50% over the ‘control’. The natural frequency, 
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Figure 6.7. The frequency response of the sartorius muscle obtained 
by random stimulation in a 'control' period (0) and 
after a 5 sec tetanus at 50/sec (@) at 16°C. 
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fy» decreases from 4 to 2.9 Hz. A slight increase is also observed in 
Ce | 

From the studies described in the previous section, it was 
observed that maximal potentiation defined by Py usually occurred in 
the first few seconds after a tetanus for temperatures greater than 
| 25°C and at 30 to 40 sec at lower temperatures (below 10°C). In the 
present experiment when a random train of 30 sec duration was applied 
after a tetanus, tension was recorded when post-tetanic potentiation 
might not be maximal, if the temperature were below 10°C. Even with 
this duration, the stimulus train produces some potentiation of the 
isometric twitch following the ‘control! random stimulation. The 
effect on the isometric twitch is however small compared with that 
produced by a tetanus. A random train of longer duration would bring 
about an early fatigue, especially at higher temperatures. 

Four muscles different from those described in Section 6.1 
were used to examine the effects of tetanus and temperature on the 
frequency response. Figure 6.8 shows the second-order parameters 
Go, fn and t of the frequency responses obtained from a typical experi- 
ment. Several results follow from Figure 6.8: 

(1) In the ‘control frequency response, raising the 

temperature decreases Gg and t, and increases fy. 

(2) Gg is greater in the 'post-tetanic' response; the 

percentage-increase becomes smaller as temperature 
is raised. 

(3) fy, is smaller in the 'post-tetanic' response; the 


percentage-decrease is about 30% at all temperatures. 


pena aca it Vo ater inke 
motighdetind aldvtst-tewg nem mente 25) eee Alle 
Abn OND” MOE iNet: Ataris ai, 9) «Hate ee vin a 
we Mis rte ee srabe ra rian aulubitee subt swath, att 
BAT ini teladirenyiIns “Levies! aah vibes Vipt feat et vrais 
Pet) APNG eit Ti hee “voenhcn zt cat! hha VTL PE oo faethe 
grade lice anearnt mgrol lo niet sours ~eainer gs ® vd one bong 
ae Ayre ig! Isls To ¥ i le fomlse weet vicse@ ea (pode: 
Lai iy) Fait oil badlinmazel ay art mi? vrei! traiaeliselly Sua 
no aura Taab bil einsysl I> or ae hy air candane oF tney One 
eiwiederiict THEM baee Sd aite Fo ocyy)s .. saeRatewy comeupett 
seadne (eo layh 4 Rant ban ePlo ead > yieulipwTh-ott ta 4 Ons gy a 
PRP saihitner) wal lal abiughy |weaved tren 


SP orl ale (setetds: “omnes! berince! eit) 11) 
raf Quetta! tn .o nr 99 2gmberset whem, 
etna ae ental re . 14) 


7 
” 


af 


a 
- 


- 


- 
Ty 


7 


: 


Figure 6.8 


73 


fa (Hz) 


0 He) 20 
Temperature (°C ) 


The influence of tetanus and temperature on the second- 
order parameters: Gog, zero-frequency gain; f,, natural 
frequency, and ~, damping ratio (shown in ordinate) used 
to describe the frequency response of the sartorius 
muscle obtained by random stimulation of the muscle nerve. 


0, 'control' and @, post-tetanic frequency response 
parameters. Data from one muscle. 
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(4) A nearly linear relationship exists between fy, 
and temperature. 

(5) Differences in t between the 'control!' and the 
"post-tetanic' response become marked at higher 
temperatures. 

These results are representative of the experiments per- 
formed. This is also true whether temperature is raised or towered. 
At higher temperatures, Go was sometimes found to be smaller in the 
'post-tetanic' response. This was probably due to muscle fatigue. 

At lower temperatures, Gp was always greater in the 'post-tetanic' 
response; ¢ did not change appreciably although both TZ. and Tir were 
increased after a tetanus. The coherence function is greater than .8 
for most of the frequencies(in the frequency response); this indicates 
the nerve-muscle system is nearly linear with this mean rate (5/sec) of 
stimulation. 

The most consistent changes were observed in f,. The 
rates of change of tension defined by 1/Tg and 1/TLip vary almost 
linearly with temperature (see Figure 6.5). Similar relationship is 
also found between fry and temperature, as is obvious from Figure 6.9. 
The data were obtained from four muscles and were fitted by a straight 
line (corr. coeff. = .98). Similar results were obtained for fy in 
the 'post-tetanic' response, although it was decreased at any temp- 


erature and the rate of change with temperature was also smaller. 
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20 30 
Temperature (°C) 


0 lO 


The relation between temperature and natural frequency, 
fy.) 01. the Sartorius muscles athe tT idsedsl ine inierceapTs 
the ordinate at .12 Hz (0°C); f}, increases at a rate of 
1.5 Hz/10°C. Data from four muscles. 
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CHAPTER 7 


DISCUSSION 

The characteristics of an isometric twitch are manifestations 
of the underlying chemical and mechanical processes. Hill (1951c) 
suggested that the maximum tension was reached as a balance between two 
opposing processes, internal shortening (mechanical) on the one hand, 
and decay of the active state (chemical) on the other. In the model 
presented here, the sliding filament theory forms the basis of the 
mechanical process while an activation mechanism depending on ayy 
constitutes the chemical process. The parameters governing the two 
processes in the model are many and are subjected to change under 
different circumstances. In this discussion, an attempt is made to 
elucidate the relative importance of the parameters from their influence 
on the characteristics of the isometric twitch, and to relate them to 
the structural and functional differences observed in skeletal muscles. 
The effects of tetanus and temperature on the characteristics of an 
isometric twitch observed in the studies here and elsewhere are also 


discussed in light of the model and other experimental evidence. 


7.1 Parameters cf the Sltding Filament Model. 

The parameters in the model have different, and characteristic, 
effects on the form of an isometric twitch. A comparison of Figures 4.3 
and 4.4 show that the constraints f,, g, and go for the making and 
breaking of cross-bridges have a smaller effect on the time course of 


the twitch than the parameters in the activation factor. f) has a 
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Similar but greater effect on the contraction phase than g,, while 
they have opposite effects on the relaxation phase. This seems 
surprising because f,; represents the making, while g, the breaking 
of cross-bridges. One possible explanation is as follows. In con- 
sidering the time course of tension development, we must take into 
account both the rate of making of cross-bridges and the number of 
cross-bridges available. A faster rate of breaking of cross-bridges 
(a larger g,) will mean that a greater number of i ea fi are 
available for making. The greater rate of cycling of cross-bridges 
permits the series elasticity to be extended more rapidly and hence 
lead to a greater rate of tension development. The relaxation phase 
is, however, a true indication of the opposite nature of ff; and gq}. 
On the other hand, go has negligible effect on the form of a twitch. 
This is because the value of go is large compared with either f, or 
Oi, SOntnal aechange by.a@ faction of foun has little overallsetieci on 
the time course of a twitch. The number of cross-bridges made in the 
region where go has a non-zero value is small, and the probability of 
the cross-bridges to be made in this particular region is also low 
under the isometric condition (Huxley, 1957). Although internal 
shortening occurs in the present model, and it is thus not a truly 
isometric condition, the velocity of shortening is stil! toosmal|! to 
have a great effect on the interaction of cross-bridges. The net 
effect is that f1;, the rate of making of cross-bridges has a greater 


influence on the characteristics of the isometric twitch than g}. 
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7.2 Parameters of the Acttvatton Factor. 

The activation mechanism is developed for frog sartorius 
muscle at a temperature of 0°C. Unless stated otherwise, any change 
would be discussed in reference to the 'standard'! activation and 
"standard' muscle at this temperature. 

Among the parameters in the activation factor, 8 has the 
greatest effect of al parameters on the relaxation phase of the 
isometric twitch. eae assumed to be the rate of release of Catt 
from the bound state in the myofilaments. I+ is generally believed 
that relaxation results from the removal of Catt from the myofi lament 
space by a process which resembles Catt-chelating by EDTA (Bozler, 
1954; Watanabe and Sleater, 1957; Ebashi, 1960, 1961b), or by a 
"relaxing factor' associated with the sarcoplasmic reticulum if the 
latter structure exists eons eee 1964). Therefore, 8 can also 
be assumed to be the rate of removal or re-uptake of cai by the 
sarcoplasmic reticulum. Whether a clear distinction can be made 
between the two processes is not known. 

Experimental evidence suggests that potentiation of the 
fsometric twitch pestis from the prolongation of the active state 
(Goffart and Ritchie, 1952; Ritchie and Wilkie, 1955). In the model, 
decreasing 8 increases the magnitude and prolongs the time course of 
the activation factor (see equation (2.12) ). The time course of the 
activation factor is similar to the active state curves studied by 
many workers (Hill and Macpherson, 1954; Macpherson and Wilkie, 1954; 
Close, 1965). This permits a greater number of cross-bridges to 
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greater peak tension is attained, but the time course of relaxation 
is also significantly prolonged. 

Similarly, a , the rate of calcium binding by the myo- 
filaments, is not necessarily distinguishable from the rate of calcium 
release by a stimulus. In the model, changing a has a relatively 
smal! effect on the characteristics of the isometric twitch. The 
time course of the activation factor in the mode! is determined by 
both a and B. Since the 'standard' value of a (65/sec) is already 
large compared with that of 8 (6/sec), changing it (up to a factor 
of four) does not alter the time course of the activation factor 


appreciably, and hence the isometric twitch. 


7.3 Acttvatton Mechantsms in Fast and Slow Muscle Ftbers. 

It has been te to classify muscle fibers according to 
their speed of contraction, and relaxation; their motor innervation; 
and in some cases to their histochemical properties. Histochemical 
studies have revealed there are in general three types of extrafusal 
fibers (fast, intermediate and slow) in mammalian striated muscles 
(Nachmias and Padykula, 1958; Stein and Padykula, 1962; Romanul, 
1964; Henneman and Olson, 1965). Studies at the motor unit level have 
shown that a smaller contraction time is associated with a larger 
diameter of the motor axon and muscle fiber (Wuerker, McPhedran and 
Henneman, 1965; Appleberg and Emonet-Denand, 1967; Burke, 1967). 
Mammalian extrafusal muscles are generally composed exclusively of 
twitch fibers, which conduct action potentials. 


Amphibian skeietal muscle is generally regarded as being 


ya Maen beet ‘oP eanelt Ae toad eit 
eee eee ‘hahaa? oar aids pen 
cnahaalaiin Sapte RA 6 rate bones aga 
Height pettalvntcn Ge twnees soi 2 ak ee en ts, Pest 
ati shereual ont SOR Bom eaten | 


a 


| : 
GT ip aye Lore wit. oy We babene gh RUS 
ar gritos ANT Biscam 9) ini gr feied oem eae) #, 7 


postage visti, Vad fiat oa A tie hot isatineg be esome tlelit > 


iagioabteet es? sen eager ee no, ; 
‘ietbatea ts Geer aah! tabditby 2 a chet SEL Aeier seat velbare  “y 7 
ee) ean Bia) vid be (idere n) sedis Gas cpdigeidtel (reek) oaedie. an » 
ee - 


ee shat . 
ies ceamapie! tee ears a Ph 


80 


composed of two types of muscle fibers: 'fast' and 'slow'; these are 
in turn innervated by the large- and small-nerve motor system (Tasaki 
and Mizutani, 1944; Kuffler and Gerard, 1947). A fast fiber can con- 
duct action potentials and produce a twitch, whereas slow fibers show 
no regenerative electrical activity and respond to local depolarization 
with a local contraction (reviews: Peachey, 1961, 1968). Slow fibers 
also have a lower fusion frequency and produce a tetanic tension which 
is better maintained than the fast fibers. Lannergren and Smith (1966) 
described five fiber types in the iliofibularis muscle of Xenopus laevis 
although the different types were classified into three main groups: 
two 'fast' and one 'slow! with overlapping functional properties. 
Studies on the motor units in the same muscle have shown both the 
‘fast! and 'slow' types can be subdivided (Smith and Lannergren, 1968). 
These authors also point out there seems to be a relationship between 
the properties of the motor unit and the diameter of its motor axon. 

Another difference between fast and slow muscle fibers 
resides in the sarcotubular system. It is generally thought that the 
slow fiber lacks a proper T-system (for review: Hess, 1967, 1970). 
Evidence from Flitney (1971) shows that frog slow fibers have an 
elaborate T-system network and sarcoplasmic reticulum. Their relative 
volumes are, however, about half as much as in a twitch fiber. 
Costantin, Podolsky and Tice (1967), using 'skinned' preparations, 
showed that the sarcoplasmic reticulum in a slow fiber is capable of 
accumulating calcium, although it appears to be less effective in 
doing so than that of a twitch fiber. 


Based on this evidence, two assumptions are made for the 
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activation mechanism in a slow muscle. The amount of calcium released 
per fiber volume by a stimulus is smaller, and so is the rate of cal- 
cium uptake by the sarcotubular system. In the model, this would be 
a reduced Co and 8. As shown in Chapter 4, this gives a twitch which 
has a smaller Py, a lower dP/dt, and a longer Tg and Tin. A tetanus 
which is produced by reducing Co and 8B also has a smaller maximum 
tension which is better maintained, a longer time to peak and a lower 
fusion frequency. These types of response are observed experimental! ly 
in some slow fibers (Smith and Lannergren, 1968; review: Hess, 1970). 
Another fundamental difference between 'fast! and 'slow' 
fibers is also found for the maximum rates at which the sliding fila- 
ments can interact. Measurement of the heat production during a 
tetanus (Floyd and Smith, 1971) and of the shortening velocity of the 
sarcomeres in ‘skinned! fibers (Costantin, Podolsky and Tice, 1967) 
and intact fibers (Peachey, quoted by Costantin et al., 1967) show 
that these processes are probably 10 to 30 times slower in the slow 
fiber. The myofibrillar ATP-ase activity, the property which ulTi- 
mately determines contractile performance (Barany, 1966), has also 
been shown, by histochamical means (Engel and Irwin, 1967), to be low 
in the slow fiber. It is, therefore, reasonable to suppose that since 
its myosin is capable only of a comparatively slow rate of activity, 
the rate at which it must be supplied with calcium need not be so 
great as in a twitch fiber. The rate at which the sliding filament 
can interact is probably equivalent to the constraint fy, the rate of 
making cross-bridges in the sliding filament model. Decreasing fy 


does give a twitch similar to that found in a muscle which has a 
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poorly-developed sarcoplasmic reticulum and again shows that fj is one 
of the more important parameters in the model. 

There is a general, though not invariable correlation between 
speed of contraction and relaxation of striated muscle fibers, and the 
volume and organization of the T-system and sarcoplasmic reticulum 
(review: Smith, 1966). In the attempt to model slow muscle, excitation- 
contraction coupling does not enter explicitly into the model. Much 
of the behavior observed experimentally in slow muscles is, however, 
obtained just by reducing Cg and 8. Ina real muscle, it is possible 
that the slowness of the activity cycle of the slow fibers reflects 
not a radically different excitation mechanism (to the extent of Gaal 
binding and removal), but a relatively more important Cane exchange 
between the surface membrane and the exterior of the fiber (Makinose 
and Hammad, cited by Hasselbach, 1964). Even in exceptional cases 
such as the asynchronous flight muscles of certain insects (Pringle, 
1965) which have a well-developed T-system but a reduced sarcoplasmic 
reticulum, some of their behavior could be obtained by modifying an 
activation mechanism which is very simitar to that in the present model 
(Julian, 1969). | 

In the extreme case where a 'threshold' is included in the 
activation factor, the effects on a twitch or tetanus obtained by 
reducing Cg and £ are exaggerated. If only one stimulus is given, 
no contraction is produced. These responses are observed in a true 
slow fiber. The tetanus obtained by eramilaring 'SSR (reducing Co 
and B by half) + threshold’ muscle (Figure 4.8) at a rate of 20/sec 
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muscle although it still declines at a relatively fast rate. Kuffler 
and Vaughan-Williams (1953) observed that with a stimulation rate of 
4/sec at 20°C the tension record of frog iliofibularis (a muscle con- 


sidered to be 'slow') is fused and maintained for many seconds. The 
decline of tetanic tension in the model results from an absence of 
replacement of the Catt store, whereas in a real muscle the Cat? store 
is continuously replenished. Winegrad (1970) reported that only 2 

sec after a short tetanus, some of the calcium released was already 
transported back to the terminal cisternae, the main release sites. 
This feature has not been included in the model, because the process 
involving the redistribution of Catt is complex as will be obvious 
later. 

Other characteristics of a muscle are demonstrated from its 
frequency-tension curve (Matthews, Rack and Westbury, 1969). This 
curve depends to a certain degree on the way in which tension is cal- 
culated. The curves for the 'standard' and the 'sSR' muscle rise very 
sharply without a sigmoid shape (Figure 4.9). This is partly because 
of the activation itself which is modelled at 0°C. At this temperature, 
the twitch/tetanus ratio is almost one and the average twitch tension 
is relatively high. Fusion frequency is low and the tetanus is almost 
completely fused at only 10/sec for the 'standard' muscle (Figure 4.7). 
The '£SR + threshold’ frequency-tension has a sigmoid shape and rep- 
resents the usual response found in real muscles (Matthews, 1959; 


Kernell, 1966). 
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7.4 Series Elastic Element. 

The series elastic element constitutes an integral part of 
the mechanical process of contraction. Hill (1951a) noted that an 
added compliance reduced the peak tension and the maximum rate of 
tension development while contraction time was increased. The changes 
observed in the twitch by changing k are similar to those observed by 
Hill (1951a) in his experiments, although the contraction time in the 
model is much more sensitive to changes in k (Figure 4.6). In face 
of the closeness of fit in Pt and dP/dt, between experimental and 
theoretical data, the reasons for the discrepancy in Te are unknown 
at present. 

Increasing k (reducing the compliance), however, affects the 
force generation process in a different way. As shown in the Fourier 
transform of the isometric twitch (Figure 4.11C), an increase in k has 
negligible effect on zc, although Gg and fy, are slightly increased. 
Recent studies on the frequency response in the cat (Mannard and Stein, 
1972), however, showed that increasing the stiffness by stretching the 
muscle increased Gg and ¢, while f, is decreased. In a real muscle, 
there is ample evidence that stretching affects the activation mechanism. 
The prolongation of the active state by stretching may result from an 
increase in the amount of activator released, a reduction in the rate 
of removal of activator, or a combination of both (Close, 1972; review: 
Sandow, 1965). In some instances, stretching a muscle can release 
sufficient calcium to produce a contraction (Armstrong, Huxley and 
Julian, 1966; Pringle, 1968; Brown, Goodwin and Matthews, 1971). The 


results from the model are as expected for a simple second-order system 
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and are also obvious from the tension curves for different k (Figure 
4.5). For example, an increased k would shorten the time course of 
contraction and relaxation; this gives a larger fy, in the frequency 
response. In the model, since the activation mechanism is independent 
of any change in k, the discrepancy of the results from those of 
Mannard and Stein (1972) seemsto lie in some kind of feedback which is 
absent in the model, but present in a real muscle. A strict comparison 
is probably not justified because the frequency response obtained by 
these authors was by random stimulation of thecat soleus muscle at 37°C 
whereas that in the model is just the Fourier transform of a twitch at 
0°C, although the second-order representation is good for either one. 
Temperature is another factor which may affect the comp!iance 
in elastic elements of the muscle. Walker (1951) found that extensi- 
bility and stress relaxation are only slightly affected by as much as 
20°C change of temperature in unstimulated rat triceps surae. These 
findings indicate that the elastic properties of resting rat muscle 
are changed very little by change of temperature. Hill (1951a) showed 
that the compliance of toad sartorius, as indicated by resistance to 
stretch, is only slightly increased by as much as 20°C decrease of 
temperature. Jewell and Wilkie (1958), in a very carefully controlled 
series of experiments, found that very large decreases of temperature 
produced only small increases of extensibility, as observed in the 
load-extension curve for series elastic elements in toad sartorius. 
It must be recognized that compliance does play a role in limiting the 
rate of tension development. However, since compliance of series 


elastic elements seems to be rather temperature-independent, it is 
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unlikely that more than a small part of the marked reduction of the 

rate of tension development induced by cooling can be attributed to 

an increase of compliance of the series elastic elements. Also, for 
changes in other characteristics of isometric twitches which follow 

temperature change, an explanation must be sought elsewhere, for 


example, in the force generation and the activation processes. 


7.5 Effeets of Tetanus and Temperature on the Charactertsttes of the 
Isometric Twttch. 
(t) Peak Twitch Tenston, Pr. 

Tetanus and temperature produce systematic changes in the 
post-tetanus twitch tension Py. The effect of temperature on the 
magnitudes of twitch and tetanic tension development of skeletal 
muscle vary widely as reported from different laboratories. Jewell 
and Wilkie (1958) observed both an increase and a decrease of twitch 
tension in frog sartorius with an increase of temperature. Winchester 
(1969) reported an increase in twitch tension from 10° to 20°C and 
a decrease above 20°C in the normal rat diaphragm muscle, while an 
increase was observed ‘from 10° to 35°C in the denervated muscle. In 
the lower range of temperatures (below°20 C), most workers report a 
consistent decrease in twitch tension, and increase in tetanic tension 
with an increase of temperature (Hill, 1951c; Kelley and Fry, 1958; 
Walker, 1949; Close and Hoh, 1968a). In the experimental studies 
reported in this thesis, it was observed that the peak tension in the 
control twitch decreases by about 49% with a change of temperature from 


3z to 30°C. However, Py in the maximally potentiated twitch decreased 
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only by 22%. This means that the tension in the potentiated twitch is 
almost doubled at higher temperatures, while at 34° the increase over 
the 'control' is only about 15%. 

Close and Hoh (1968a) observed that with rat extensor digi- 
torum longus (EDL, considered to be a 'fast' muscle), Peri netbe weoptro |!” 
isometric twitch decreased by about 40%, while that in the maximally 
potentiated twitch remained almost at the same value with a temperature 
change from 20° to 35°C. Py in the soleus (SOL, a 'slow' muscle), 
however, increased by only a few percent with the same temperature 
change. There is little or no post-tetanic potentiation in EDL at 
20°C and Py at that temperature is about the same as the peak tension 
of the maximally potentiated twitch at other temperatures between 20° 
and 35°C. The increase in tetanic tension, however, is considerable 
(about 30%) in EDL, but is very smal! in SOL with the same temperature 
increase. A similar degree of maximum post-tetanic potentiation has 
been observed in cat fast muscles at 37°C and frog twitch muscle at 
about 20°C and has been shown to occur in single muscle fibers (Brown 
and von Euler, 1938; Ramsey and Street, 1941). A possible explanation 
suggested by Close and Hoh (1968a) for the differences in temperature 
dependence of frog twitch muscle and rat EDL muscle on one hand and 
rat SOL on the other is that a decrease in temperature decreases the 
intrinsic speed of shortening and the rate of removal of activator in 
all these muscles, thereby causing an inversely proportional increase 
in Py (Close, 1965). The decrease in temperature, however, also 
increases P+ of frog twitch muscle (for example, the sartorius) and 


adult rat EDL, in much the same way that repetitive stimulation in- 
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creases it, by increasing the degree of activation of individual muscle 
fibers with little or no change in the time course of the response. 

The above scheme offers a good explanation for the results 
obtained here and elsewhere. The time course of the maximally poten- 
tiated twitch is, however, significantly prolonged by a tetanus at any 
temperature in the range 34 to 30°C in the present experimental studies 
(Figure 6.5). In the scheme of Close and Hoh (1968a), there would 
probably be little or no change in Tg, Tin and dP/dt. The maximum rate 
of tension development was not measured, and it is doubtful whether 
the discrepancy is a result of stimulation of the muscle nerve, and 
not the muscle directly. Connolly, Gough and Winegrad (1971), by 
stimulating the sartorius muscle of the frog directly, observed an 
increase in both 7, and dP/dt in the maximally potentiated twitch. The 
meaning of the term 'degree of activation’ is rather flexible; ina 
muscle it would be proportional to the number of active fibers (Close 
and Hoh, 1968a) whereas in a single fiber it would probably be a 
measure of the capacities to bear a load and to shorten. The 'degree 
of activation' may be related to the intrinsic speed of shortening 
in a particular muscle or muscle fiber at a particular temperature. 

In light of the experimental evidence presented here or elsewhere 

(for example, Connolly, Gough and Winegrad, 1971), it is certain that 
the change in the time course of a twitch after a tetanus can be ex- 
sfathet in some other way. It does not conflict with Close's (1965, 
1968a) suggestion that the intrinsic speed of shortening is not changed 
by a tetanus and is also independent of the 'degree of activation’. 


In relation to the model, it is obvious that decreasing 8 
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gives a twitch which resembles the maximally potentiated twitch 
observed experimentally (Figure 4.1). The significance of 8, the 
rate of re-uptake of Cale is discussed in more detail later in 
reference to changes of the relaxation phase of an isometric twitch 
after a tetanus. It ts difficult to speculate about the effect of 
temperature change on the parameters of the model. Reducing 8 alone 
does give a twitch whose Tin is greatly prolonged, but T, is not 
significantly increased, and there is a very small increase in dP/dt 
(Figure 4.1), whereas the time course of the twitch of a real muscle 
is greatly prolonged by a decrease in temperature. The other rate 
constant, a, which also governs the time course of the activation factor, 
has a relatively small effect on the twitch. 

The rate of making of cross-bridges, in a loose sense, can 
be related to the intrinsic speed of shortening, which is decreased 
when temperature is lowered. From a look at the separate effects of 
8 and f,; (Figures 4.3 and 4.4 respectively), it is expected that 
reducing both 8 and f; would give a response quite similar to that 
of a real muscle with a decrease in temperature. However, temperature 
affects all the parameters in a complicated manner and experimental 
studies seem to be the most rewarding approach in this regard. 

In studying the effects of tetanus on the characteristics 
of the isometric twitch, Connolly, Gough and Winegrad (1971) observed 
that the potentiation was not maximal until some time after the tetanus. 
He suggested that this is due to the net effect of two processes. 

The fast process is responsible for the delayed increase in Pz 


and a greater 7, after the tetanus. This is aiso identified with a 


7 


as _ 
ine 4 ut? .9)9eag a gale rt. 
ni i) bem ca warts ioe a tat 


» | 


wae a 7 wa ; 
Pail. ore curt (tb were Eee 
7 ; foe a 
_ hee ee kn 7 - 
~ 4 . eo / = in 7 14 7 2 af ’ " ‘ 
okial etek aad 1 x) + e « ye . + y _ 
| »¢ >| we re | Ti - ae ty 46 Ce wt “a ( i. " ras ” wt ess ~ fon 
- 7 7 airs ’ 
i ty Piads at f ts A 214) oe nan 4 1) oot eo %t! 7 
> : - ' Ss 
“te 4 } Cc. Br ’ 2tg2greanec *# onpey i> Oh 3 wouter . =, 
y ; rw 0 : 
a > 7. S 
! \ "hh be : “Wes aig Oh, dd. Richie sus on, 
= i . le 
a a 
: aa = ' y! tans h iiegle 
: ' e 
: _ ' 
| ¢ j ~ ~~ i : = he) aper¢ fiw tee egy 
: we 
s » 2 ‘ o 
= * ‘ 
¥ re ' 1? ¢ i 
7 . , »* 7 ip wt 
07 ’ ¥ ee | 
4 ‘i aed : i 
17 ° (| : 1 a 
iyi seit rai? nay ) ‘ 


AUP ah petre Gartiaias |! finu rer! var “in bw PE 
Merit’ 4 | is. 
oy al 1% Ow’ ts i qe ei ay) 7 “ft & b's os 
- ana 
tty wat 9 J " le “ 


a AF ts vps fen 
se, 


' ‘ > 
- , LP Oe ; 
ar I as Oo: oes 


_ 


90 


fast-diminishing inhibition: during the tetanus, the Ca*t stored in 
the terminal cisternae is depleted. As a result, the amount of cal- 
cium released immediately after a tetanus by a stimulus is smaller. 
This would result in a smaller Fee (nts tnhibttion ts -alileviared 
when the calcium store is replenished by relocation of Catt from the 
longitudinal tubules and intermediate cisternae portions of the 
sarcoplasmic reticulum after re-uptake. These authors analysed 

the change in Pz, dP/dt and Tg and found that these characteristics 
followed a time course and had a Q)9 similar to that observed for 
the calcium movement in the sarcoplasmic reticulum after a tetanus. 
The fast-diminishing inhibition results from the depletion of 
calcium and from its fast replenishment. 

In the studies reported here, by stimulating the muscle 
nerve, similar changes were observed for Pz. In the temperature 
range 33° to 25°C,Pt usually increased nearly linearly to a maximum 
and then started to decay exponentially. This is expected if Pz 
indicates the sum of a fast-decaying and a slow-decaying process. 

At higher temperatures, the potentiation reached a maximum in the 

first few seconds after a tetanus. Analysis showed that the rising 
phase has a @19 of 2./ compared to 1.4 as found by Connolly, Gough 
and Winegrad (1971). I+ takes about 15 sec to reach its peak at 20°C. 

The falling phase of the potentiation is associated with a 
slowly-decaying potentiation process immediately following a tetanus. 
These authors, however, reported that the underlying processes res- 
ponsible were very complex but did not elaborate on their exact 


nature. In the studies reported here, the slow phase as measured from 
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the falling phase of potentiation (using Pz) has a half-time of about 
1z min at 20°C and a Qi9 of 2.9. 

The amount of calcium released from the terminal cisternae 
is large during the tetanus. This large amount would tend to saturate 
the process of re-uptake by the sarcoplasmic reticulum. It is then 
not unreasonable to suggest that this slowly-decaying potentiation is 
due to a decrease in the rate of re-uptake (8) of Catt to a minimum 
value after the tetanus. The fact that Tin always has the greatest 
value immediately after a tetanus supports this view. 

The changes in P;, and Tiy after a tetanus can also be 
associated with some of the parameters of the model. If Co is 
reduced, the peak tension is also reduced as observed in the experi- 
ments (Figure 4.1). Reducing 8 in the model prolongs the relaxation 
phase but. also increases Pt (Figure 4.1). Pz is, however, mostly 
dependent on Cg whereas the relaxation phase is governed mainly by 
B. At higher temperatures, the Can stcrents continuously and 
quickly replenished, even during the course of the tetanus. As a 
result, Cg may not change significantly after the tetanus; whereas 
since the return of 8 to its pretetanic value is a much slower 
process, its effect is exaggerated, leading to a maximal potentiation 
almost immediately after the tetanus. The decay of this potentiation 
would be a total replenishment of the calcium store, with 8 increasing 


from its post-tetanus value. This process would take more than 15 


minutes at lower temperatures. 
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(tt) Contraction Time, Te and Half-Relaxation Time, Thy. 

The reciprocals of Tg and Ti» are referred to as the 
contraction-rate and relaxation-rate respectively. Both the con- 
traction and relaxation rates are observed to be minimal after the 
tetanus but the time course of recovery is more variable than Py. 
From the slopes obtained by fitting the points with straight tines 
in Figure 6.5, the ratio contraction-rate/relaxation-rate (equal 
to T!y,p/Tg) equals 1.5 and 2 for the 'control' and the maximal ly 
potentiated twitch. The ratio Tir/Te is, in fact, greater for the 
potentiated twitch at higher temperatures, whereas at lower temper- 
atures, this ratio does not change appreciably. Supporting evidence 
was also obtained from the frequency response measured after a 


tetanus, as will be obvious in the discussion later. 


7.6 Effects of Tetanus and Temperatures on the Frequency Response 
of the Isometric Muscle. 

The representation of the muscle in these studies by a 
second-order system is extremely good, as shown by analysis in the 
frequency domain. This is true whether the Fourier transform was 
computed from a twitch (Figure 6.6) or when the frequency response 
was obtained by random stimulation (Figure 6.7). The coherence 
obtained in the latter analysis is high (> .8 for most frequencies). 
This shows that with the low rate of stimulation used (mean rate = 
5/sec), the muscle responds quite linearly. At higher mean rates, 
the coherence becomes more variable (Stein, French, Mannard and Yemm, 


1972)% 
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The percentage increase in zero-frequency gain Go is 
greater in the frequency response obtained after a short tetanus 
of a few seconds duration at lower temperatures; here the fusion 
frequency is lower while the muscle does not fatigue quickly. As 
a result, potentiation lasts longer and gives a larger Gop in the 
frequency response. At higher temperatures, Gg is more variable. 
Sometimes a decrease is observed although a small increase is most 
common. This seems paradoxical because potentiation occurs at all 
temperatures and should give a larger Go. In the studies on the 
effect of tetanus on the characteristics of the isometric twitch, 
it was, however, observed that Pemay be depressed atterm a poten= 
tiation (compared with the pretetanic value). The depression of Pz 
seems to be related to the smaller Gg obtained sometimes. This may 
result from muscle fatigue or from some other unknown causes. 

The most consistent changes with temperature are observed 
in fye This is partly because Tg and Tip» always increased after a 
tetanus, although P+ sometimes decreased to a value smaller than the 
‘control’ after a potentiation. In fact, a nearly linear relationship 
(cortemcoctt ws —..90) lis found ifog foainiane || control! eel nihe.Ussec 
after a tetanus, f,, is decreased for every temperature but also varies 
nearly linearly with temperature. 

Hutchinson, Koles and Smith (1970) found that in toad mye- 
linated nerve fibers the conduction velocity varies linearly with 
temperature. If f,, indicates the speed of response of the muscle, 
it seems that at least one neuronal and another mechanica! response 


in the neuromuscular system vary linearly with temperature. Whether 
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this linearity is important to movement co-ordination in poikilothermous 
animals at different temperatures is at present unknown. 

At lower temperatures, ~ was usually somewhat greater than 
1, but a tetanus did not seem to change its value appreciably. At 
higher temperatures, the effect of tetanus was very marked and &% was 
always increased. Since at any temperature a tetanus increased both 
‘ies and Tins the above findings seem paradoxical. This can, however, 
be explained as follows. If the impulse response is calculated for a 
second-order mode! with different values of z, the measured ratio 
Tir/Te can be related to x. For example, it is found that T1,/T, 
= 1.9 and 1.68 for c = 1.2 and 1 respectively. In other words, z 
depends on the ratio T!y/Tg, and not their absolute values. At all 
temperatures, T,, Tiy and the fanro Tiy/Te are increased by a 
tetanus, but at higher temperatures, it was observed that a tetanus 
brought about’ a disproportional increase in fl, relative to To, 
and therefore ¢ is increased markedly. Because f,, Is observed to 
increase with a temperature increase, in both the 'control' and the 
"post-tetanus' response, it is in a sense increasing in proportion to 
both 71, and Tg, and is thus a good indicator of the "sluggishness'! 


of the response of the muscle. 


In conclusion, the second-order representation of the 
muscle and the mode! studied here has been demonstrated to be extremely 
good, as was also reported in mammalian muscle (Stein, French, Mannard 
and Yemm, 1972; Mannard and Stein, 1972). These authors suggested 


that, despite the many rate constants governing the chemical and 
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mechanical processes underlying activation and contraction, only two 
are rate-limiting. In the work reported here, it was also found that 
among all the parameters in the model, two seem to affect the char- 
acteristics of the isometric twitch most significantly. These are 
fi, the rate of making of cross-bridges, and 8, the rate of re-uptake 
of calcium by the sarcoplasmic reticulum. These rates also seem to 
determine and can be used to explain to some extent the functional and 
structural differences found in real muscles and the changes of the 
tension response under different physiological conditions, such as 
after a tetanus and temperature change. Additional rate constants 
seem to determine the changes in B after a tetanus and the movement 
of Catt in the sarcoplasmic reticulum, and between the latter 
structure and the myofilament space. These can only be elucidated 

by studying the biochemical and biophysical reactions underlying 

the release, binding and re-accumulation of calcium. Other reactions 
involving the utilization of ATP, which is closely related to the 
chemical process of calcium activation can be included in the modei. 
The difference in activation-contraction coupling mechanisms observed 


for fast and slow muscles could also be included in future work. 
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